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ABSTRACT 



This report details the development and provides the 
documentation for custom computer software applications that 
evaluate antenna parameters. The applications are written in 
Mathcad 3.1 for the following antenna types: linear, planar, 
and circular arrays; folded dipole; caged dipole; parabolic 
reflectors with helical and spiral feeds. Inputs to the 
applications are limited to the antenna's physical dimensions. 
In some cases, ground parameters are required. 

The chapters are structured to provide the user with the 
necessary information needed to use and interpret the software 
for each antenna type. The software applications are provided 
as appendices and give examples of each antenna type. 

Outputs of the applications provide various numerical and 
performance predictions, as well as far-field radiation 
patterns. The results computed are consistent with 
predictions provided in applicable publications, as well as 
those calculated by numerical antenna analysis programs such 
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as ELNEC. 
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I. BACKGROUND AND PURPOSE 



This thesis details the development of user friendly 
mathematical applications capable of computing the radiation 
pattern and other pertinent antenna parameters of an antenna 
or antenna system based on available information, primarily 
dimensional information obtained from photographs. The 
applications developed are compatible with any IBM personnel 
computer using MS-DOS version 3.2 or higher and a math 
coprocessor. The applications are written in Mathcad version 
3.1 engineering software and can be run using either Mathcad 
3.1 or Mathcad version 4.0. These applications are intended 
for use by the Naval Maritime Intelligence Center 
( NAVMARI NTCEN ) . 

As already mentioned, user inputs are limited to antenna 
dimensions based on physical measurements primarily obtained 
from photographs, although other source data can be used. The 
applications are flexible to the extent that other parameters 
may be estimated and used as inputs to increase the accuracy 
of the computations to gain better insight into the antenna's 
performance. The Mathcad applications provide various 
performance predictions as well as a graphical representation 
of the antenna's far-field radiation pattern. The necessary 
background information needed to interpret the application's 
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formulas and displays are provided in the corresponding thesis 
chapters . 

Dietrich [Ref. 1] and Gerry [Ref. 2] have completed the 
first and second reports of this project, respectively. This 
thesis is the last in a series of three reports intended to 
fulfill the NAVMARINTCEN statement of work. 
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II. INTRODUCTION 



The increasingly sophisticated design by foreign 
countries of antennas warrants the development of improved 
analytical methods to provide timely and accurate technical 
assessment of these antennas. Through the process of reverse 
engineering, foreign communication, navigation. 
Identification, Friend or Foe (IFF) , and radar antennas are 
evaluated to determine their capabilities, limitations, and 
vulnerabilities. Having limited information available 
concerning the parameters of an antenna, usually restricted to 
photographs that provide only the type and physical dimensions 
of the antenna, current methods of antenna evaluation are slow 
and tedious. 

With the availability of high-speed personal computers 
matched with sophisticated off-the-shelf engineering software, 
the assessment of antennas using only photographic information 
may now be obtained rapidly. The objective of this report is 
to provide NAVMARINTCEN with user friendly Mathcad software to 
achieve this type of performance in their evaluation of 
antenna systems. 

As with any computer analysis of antennas, the accuracy 
of the results depend upon the sophistication of software 
used, the complexity of calculations, the reliability of input 
data obtained from photographs, and the accuracy of any 
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required input data that must be estimated. This report and 
the Mathcad applications are written with existing engineering 
equations for the type of antenna analyzed and inputs are 
limited to physical dimensions. Additional information such 
as ground effects and characteristic impedance may be included 
to increase the accuracy of the analysis. To this end, this 
report should provide an initial understanding as to the 
capabilities of the antenna system in question. 

The chapters of this report document the equations and 
assumptions used for each specific antenna type and is written 
as a comprehensive reference for the software. Copies of each 
application are included as appendices to provide the user 
with a printed illustration of the software. 
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III. THE UNIFORM LINEAR ARRAY ANTENNA 



Arrays are highly directional antenna system with narrow, 
steerable beams and low side lobes used for long distance 
communication and radar systems. The linear array is one 
geometrical configuration that consists of radiating elements 
lying along a straight line. As a result of the dependance of 
inter-element spacing (d) to the wavelength (X) , linear arrays 
are inherently narrowband antenna systems. 

In this Mathcad application uniformly spaced, equally 
excited, isotropic point sources aligned on the z-axis are 
used as the radiating elements. In using isotropic point 
sources, polarization of the array is not calculated and an 
assumption is made that the antenna input resistance (R jn ) is 
equal to the radiation resistance (R r ) . This assumption 
results in calculations for gain (G) that are idealized 
values. The following inputs are needed to implement the 
software : 

N = number of isotropic radiating elements 

d = inter-element spacing between adjacent elements 

0 = direction of main lobe 

O 

f = frequency of interest 

1 = antenna feed current 
0 

Z 0 = characteristic feed impedance 
The first two inputs are obtained from visual data such as 
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photographs, and the indicates input parameters that are 
either known or are estimated. The linear array geometry is 
illustrated in Figure 3.1. 



z 




X 



FIGURE 3.1 Geometry of N-element linear array of 
isotropic point sources 
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The total electric field of a linear array in the far- 
field, neglecting mutual coupling between adjacent elements, 
is determined by the vector addition of the fields of 
individual radiating elements. It can be shown that the total 
field is equivalent to the product of the field of a single 
element, selected at a reference point, and a factor which is 
referred to as the array factor (AF) : 

E (total) = [E (element) ] x [AF] (3.1) 

In (3.1), E is the vector electric field intensity, while AF 
is a scalar quantity. This result is referred to as pattern 
multiplication for arrays of identical elements. The array 
factor is a function of the geometry of the array and the 
element excitation amplitude and phase. For uniform 
excitation, by varying the inter-element spacing (d) and/or 
the phase (/?) between the elements, the characteristics of the 
array factor, and consequently the total electric field, of 
the array can be controlled. [Ref. 3: pp. 204-207] 

This Mathcad application is written for an N-element 
linear array with equally spaced, uniformly excited isotropic 
radiating elements. In addition, element phase is assumed to 
vary linearly along the length of the array. The radiation 
pattern of a specific element type is neglected in (3.1) since 
in normal usage it will have little effect when the array 
consists of a large number of elements. However, the user of 
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this software should be able to use information regarding the 
type of element and the orientation of the elements to narrow 
the range of possible 0 Q 's. The normalized array factor (AF n ) 
is [Ref. 3: pp. 212-214]: 



AF n (Q) 



1 


(0) j 


1 


N 


sin(-| i[> (0) j 


l_ 



( dimensionless ) 



(3.2) 



where 



q (0) 


= kdcos(d) + 


P ( radians ) 


(3.3) 


P 


= -kdcos (0 O ) 


(radians) 


(3.4) 



In (3.4), 0 o is the direction of the maximum array factor 
value (i.e., for a broadside array 0 O = n/2, for an end-fire 
array with the main beam directed at 180°, 0 o = n , and for a 
phased array with the main beam directed to ±60°, 0 Q = n/3) , 
(3 is the phase excitation difference between the elements, k 
is the wavenumber (2n/X ) , and d is the inter-element spacing. 

With the linear array aligned along the z-axis, the array 
factor is a function of 0 but not of <p. The visible region of 
the array in terms of 0, i|/(0), and the inter-element spacing 
(d/A.) is: 

0 < 0 < u (radians) (3*5) 
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P-icd < *|i (0) < P+^rd ( radians ) 



(3.6) 



d/A determines how much of the array factor appears in the 
visible region as defined by (3.5) and (3.6). Although d/A 
provides guidance in determining the frequency to select for 
this Mathcad application, it should be noted that most arrays 
are designed with the inter-element spacing less than one 
wavelength, usually close to a half -wavelength. Thus, linear 
arrays are characteristic narrowband antenna systems. [Ref. 
4: p. 123] 

An observation point is considered to be in the far-field 
if all of the following are satisfied [Ref. 4: pp. 24-25]: 



The array length, Nd, in (3.8) and (3.9) is the maximum 
dimension of the antenna (the array length is assumed to 
include a distance of d/2 beyond each end element [Ref. 5: 
p.55]). The minimum distance to the far-field is found by 
selecting the largest value of (3.7), (3.8), and (3.9). 

Radiation intensity (U(0)) is related to the power 
radiated in a given direction and is independent of the 



r a 1 . 6 A ( meters ) 



(3.7) 



r k 5 (Nd) (meters) 



(3.8) 




(3.9) 



9 



distance to the observation point [Ref. 3: p. 27; Ref. 4; p. 
33]. For linear arrays U(0) is given by [Ref. 3: pp. 229- 
233] : 

<7(0) = ( AF n ) 2 (W/ solid angle) (3.10) 

Radiated power (P rad ) is the total power radiated by the 
antenna and is obtained by integrating U(0) over a surface 
surrounding the antenna [Ref. 4: p. 33]. Since the geometry 
of the linear array is aligned on the z-axis and the array, 
factor is independent of <p , P pad is thus [Ref. 3: pp. 229-233]: 

P iad = 2tc f n U( 0 ) sin ( 0 ) dd (W) (3.11) 

Jo 

Directivity (D 0 ) is the maximum value of directive gain 
where directive gain is the ratio of U(0) in a specific 
direction to the average radiation intensity [Ref. 4: pp. 34- 
36]. Hence, D q is given by [Ref. 3: pp. 229-233]: 

D 0 = - U ^ max (.dimensionless) (3.12) 

^rad 

In (3.12), U equals unity and occurs at 0 Q for the uniform 
linear array. 

Effective isotropic radiated power (EIRP) is a commonly 
used communications term that is defined as the product of 
antenna gain and total power accepted by the antenna from the 
transmitter. EIRP is determined as [Ref. 4: p. 62]: 
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EIRP = P Iad Do 



(3.13) 



(W) 



Assuming that the power delivered to the array consisting 
of N-isotropic point sources is equivalent to the radiated 
power, we can estimate the radiation resistance (R r ) of the 
antenna as [Ref. 3: p. 55]: 



Using the same logic, we obtain the input resistance (R in ) of 
the array made up of isotropic point sources as equivalent to 
the radiation resistance: 



The exact gain (G) of the array cannot be calculated with 
the data assumed, but an idealized value can be computed based 
on the assumption made to determine R jn and the assumptions 
made to calculate total antenna efficiency (e t ) . The total 
antenna efficiency is given by [Ref. 3: pp. 44-45]: 



R 



I 




(3.14) 



R in = R r (Q) 



(3.15) 



e t = e r e cd = 1 — | r | 2 (dimensionless) 



(3.16) 



where 



t cd = 1 (dimensionless) 



(3.17) 



€ r = 1 — | r | 2 (dimensionless) 



(3.18) 
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( dimensi onl ess) 



(3.19) 




O 



O 



In (3.17), e cd is the conduction and dielectric efficiency 
which is unity for an ideal lossless antenna. In (3.19), Z o 
is the characteristic impedance of the transmission line and 
r is voltage reflection coefficient. Therefore, the total 
antenna efficiency is equal to the mismatch efficiency of 
(3.18) and the ideal gain for a linear array is [Ref. 3: pp. 
43-44] : 



As a result of constructing the linear array with 
isotropic point sources, the polarization loss factor (PLF) 
and antenna polarization are not calculated. Therefore, the 
maximum effective aperture (A em ) of the array when used as a 
receiver is calculated assuming negligible polarization 
mismatches. The A em is then [Ref. 3: p. 63]: 



G = t t D 0 ( dimensionless ) 



(3.20) 



G(dB) = 10log 10 (c c D a ) ( dB ) 



(3.21) 




(3.22) 



where 



PLF = 1 ( dimensionless ) 



(3.23) 
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The assumptions used in finding R r and A em result in an 
idealized value in calculating the effective height (h^) . The 
effective height is [Ref. 6: p. 42]:] 



Km = 2 . 



R i A em 



\ Tlo 



[meters) 



(3.24) 



where r} o is the characteristic impedance of free space (120tt). 

Two displays are produced in this Mathcad application. 
The first display is the rectangular representation of the 
magnitude of the array factor from -n to n . The second 
display is the polar representation of the magnitude of the 
array factor. 

The Hansen-Woodyard end-fire array was not addressed in 
this application nor were arrays of uniform spacing but 
nonuniform amplitude (i.e., Dolph-Tschebyschef f array). In 
addition, nonuniformly spaced, uniformly excited arrays were 
not addressed by this application. 

The results obtained with this application are identical 
to those found in [Ref. 3: pp. 216-240] with respect to array 
factor patterns and directivity. Table 3.1 is a comparison of 
data from [Ref. 3] for a 10-element equally excited, uniformly 
spaced, end-fire linear array while Figure 3.2 is a comparison 
of the array factor pattern of the same array. Of note, the 
array factor patterns from this application and the patterns 
from [Ref. 3] are identical as expected. 
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Table 3.1 End-fire Array Data Comparison 



ANTENNA 


REFERENCE 


CALCULATED 


DIMENSIONS 


DIRECTIVITY (D 0 ) 


DIRECTIVITY ( D Q ) 


N = 10 






d = A/4 m 


10.0 dB 


10.003 dB 


o 

o 

II 

o 

<D 






'CO 

II 

1 

a 
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FIGURE 3.2 Array factor pattern of 10-element equally 
excited, uniformly spaced end-fire linear array with d = 
1/4 m, 1 0 = 0°, and /3 - -kd 
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Table 3.2 is a comparison of data from [Ref. 3] for a 10- 
element equally excited, uniformly spaced, broadside linear 
array, and Figure 3.3 is a comparison of the array factor 
pattern of the same array. The array factor patterns from 
this application and the patterns from [Ref. 3] are identical. 



Table 3.2 Broadside Array Data Comparison 
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-1 - 0.8 - 0.5 - 0.4 - 0.2 0 0.2 0.4 0.6 0.8 1 



FIGURE 3.3 Array factor pattern of 10-element equally 
excited, uniformly spaced broadside linear array with d = 
1/4 m, 0 O = 7r/2, and (3 = 0 



17 



IV. THE PLANAR ARRAY ANTENNA 



Planar arrays, frequently used long distance 

communications and radar systems, exhibit characteristics 
analogous to those of linear arrays but have additional 
capability to control and shape the radiation pattern. Planar 
arrays are more versatile than linear arrays, providing 
symmetrical patterns with lower side lobes, and can be used to 
scan the main beam of the antenna array toward any point in 
space [Ref. 3: p. 261]. As with linear arrays, planar arrays 
are narrowband antenna systems due to the dependance of the 
inter-element spacing (d„ and dj on the wavelength (A.) . The 
following inputs are required analyze a planar array: 

M = number of isotropic radiating elements in the 
x-direction 

N = number of isotropic radiating elements in the 
y-direction 

d x = inter-element spacing between adjacent elements in 
the x-direction 

d y = inter-element spacing between adjacent elements in 
the y-direction 
f = frequency of interest 
0. = direction of main lobe at 0 = 0 * 

O 0 

0 O = direction of main lobe at <p = <p 0 
I 0 = antenna feed current 



18 



Z Q = characteristic feed impedance 
The first four inputs are data obtained primarily from visual 
data such as photographs, and the * indicates input parameters 
that are either known or estimated. The planar array geometry 
is illustrated in Figure 4.1. 




FIGURE 4.1 Geometry of M x N element planar array 
of isotropic point sources 
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This Mathcad application is written for an M x N planar 
array with uniformly spaced, equally excited, identical 
radiating elements positioned on the x-y plane. This 
application may also be used to approximate the array factor 
and field patterns for arrays that have elliptical or circular 
geometries by assuming an ellipse is equivalent to a rectangle 
or a circle is equivalent to a square, respectively [Ref. 6: 
p.187]. Since isotropic point sources are assumed as array 
elements, the polarization of the array is not calculated, and 
an assumption is made that the antenna input resistance (R jn ) 
is equivalent to the radiation resistance (R r ) . This 
assumption results in values for gain (G) that are idealized. 
As with the linear array, the radiation pattern of a specific 
element type is neglected in the analysis of the planar array 
since in normal usage it will have little effect for a planar 
array with a large number of elements. However, as with the 
linear array, the user of this software should be able to use 
information regarding the type of element and the orientation 
of the elements to narrow the range of the possible 0 o ,s and 

V s - 

The normalized array factor (AF n ) is calculated for far- 
field observations neglecting mutual coupling between adjacent 
elements with equal amplitude excitation and is [Ref. 3: pp. 
260-263 ] : 
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AFje,<t>) 



(4.1) 



sin 


(I'M 0 .*)) 




sinj 


ty( 0 /4>) ) 




M sin 


' ^ x ( 0 '4>) ' 
2 


j| 


„ . [ 4fy( 0 /4>) ' 

N sxnl — - — ^ 


) 



where 

ilr x (0,<t>) = kd x sin ( 0 ) cos ( <(> ) + P x (radians) (4.2) 

ilr y (0,<|>) = kd y sin (0 ) sin (<{>) + P y ( radians ) (4.3) 

P x = -kd x sin ( 0 O ) cos ( <J> 0 ) ( radians ) (4.4) 

P y = -ied y sin(0 o ) sin(4> 0 ) ( radians ) (4.5) 

Equation (4.1) is the normalized array factor that is a 
function of both 0 and 0 where (4.2) and (4.3) are the array 
factor phase shifts in the x and y directions, respectively. 
(4.4) and (4.5) are the progressive phase shift between 
adjacent elements in the x and y directions, respectively, 
with the main beam directed along 0 = 0 O and 0 = <p Q . In (4.2) 
through (4.5), k is the wavenumber (2n/X) . 

When the inter-element spacing between the elements (d x 
and d y ) equals or is greater than the wavelength (X) , multiple 
maxima of equal magnitude are formed. To avoid creating 

multiple maxima, the inter-element spacing must be less than 
a wavelength at the frequency of operation [Ref. 3: p. 263]. 
With the dependance of inter-element spacing on the 
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wavelength, planar arrays are characteristic narrowband 
antenna systems. 

An observation point is considered to be in the far-field 
if all of the following are satisfied [Ref. 4: pp. 24-25]: 

r £ 1 . 6 X (meters) (4.6) 



r 2 5 J ( A Tcf x ) 2 + ( ~Nd ~ y ) 2 ( meters ) 



2 ( ( Md ) 2 + (Nd y ) 2 ) , . 

r ^ i - z - ( meters ) 

A 



(4.7) 



(4.8) 



The length of an array is assumed to include a distance d/2 
beyond each end element [Ref. 5: p. 55]. Therefore, the 
maximum dimension of the planar array is: 

^ ( Md x ) 2 + ( Nd y ) 2 ( meters ) (4.9) 



The minimum distance to the far-field is found selecting the 
largest value of (4.6), (4.7), and (4.8). 

Directivity (D o ) for planar arrays is determined 
approximately by [Ref. 7: p. 78]: 

AF a AF* n L x 



D n = 4n 

° ^ 2k rn/2 



£ ” f* (AF n AF„) sin(0) dQd<\> 



(W 



(4.10) 



The maximum value of the numerator in (4.10) is unity at 8 = 
0 O and 0 = 0 o since the array factor is normalized. The 
denominator is the radiated power (P rad ) of the antenna where 
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P rad obtained by integrating over all angles around the 

antenna above the x-y plane (0 < 0 < n/2) [Ref. 3: p. 28; Ref. 
7: p. 78]. This assumes the array is always pointed toward 
one side or the other of the x-y plane. 

Effective isotropic radiated power (EIRP) for planar 

arrays is defined as the product of antenna gain and the total 
power radiated. EIRP is determined as [Ref. 4: p. 62]: 

EIRP = P rad D a (W) (4.11) 

Assuming that the power delivered to planar arrays 
constructed of M x N isotropic radiating elements is 
equivalent to the radiated power, we estimate the radiation 
resistance (R r ) of the array as [Ref. 3: p. 55]: 

= 2 j - f r | a / ) (Q) (4.12) 

I -To I 

Since the radiating elements are assumed to be ideal point 
sources, the input resistance of a planar array (R in ) is 
assumed to be equivalent to the radiation resistance: 

R in = R x (Q) (4.13) 

The gain (G) of a planar array cannot be calculated given 
the information assumed to be available, but an idealized 
value can be derived. Another factor affecting the gain is 
the assumption made when determining the total antenna 
efficiency (e t ) . The total antenna efficiency for an antenna 
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is [Ref. 3: pp. 44-45]: 



e t = e r e cd = l-|r | 2 (dimensionless) 



(4.14) 



where 



t cd = 1 ( dimensionless ) 



(4.15) 



e r = 1 - | r | 2 (dimensionless) 



(4.16) 



p = 3m — £2 (dimensionless) 

R i r, + Z n 



(4.17) 



In (4.14), e cd is the conduction and dielectric efficiency, 
which is unity for an ideal lossless antenna. In (4.17), Z Q 
is the characteristic impedance of the transmission line and 
T is the voltage reflection coefficient. Therefore, the total 
antenna efficiency is equal to the mismatch efficiency of 
(4.15) when ideal lossless antenna is assumed. The ideal gain 
of a planar array in dimensionless and decibel quantities is 
then [Ref. 3: pp. 43-44]: 



The maximum effective aperture (A em ) for a planar array 
is derived by assuming polarization mismatches are negligible 
(this application is constructed using isotropic radiating 



G = t t D 0 (dimensionless) 



(4.18) 



G(dB) = 10 log 10 (e t D o ) (dB) 



(4.19) 
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elements for which polarization loss factor (PLF) and antenna 
polarization are not calculated) [Ref. 3: p. 63]: 



A 



'em 




(4.20) 
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PLF = 1 ( dimensionless ) 



(4.21) 



The effective height (h em ) of the planar array is now 
estimated by [Ref. 6: p. 42]: 



In (4.22), rj 0 is the characteristic impedance of free space 
(1207T) . 

Four two-dimensional array pattern displays are produced 
in this Mathcad application: the first is in the x-z plane; 
the second is in the y-z plane; the third display generated is 
a plane perpendicular to the x-y plane when 0 = 7r/4 ; the last 
display is a plane perpendicular to the x-y plane when <p = <p 0 . 
The last display is redundant when <p is either 0, 7r/4, or n/2. 

Results obtained with this application are the same as 
those found in [Ref. 3: pp. 260-274]. Directivity was 

calculated using the approximate equation (4.10). An 
alternative approximation for directivity is given by either 
[Ref. 3: p. 272-273]: 



h 




(4.22) 
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( dimensionless ) 



(4.23) 



D 0 = n cos ( 0 O ) D x D y 



or 

2 

L> “ ( dimensionless ) (4.24) 

Q a ( rads 2 ) 

In (4.23), D x and D y are the directivities of broadside linear 
arrays of length and number of elements Md x , M, and Nd y , N in 
the x and y directions, respectively. In (4.24), fl A is the 
beam solid angle. 

Table 4.1 is a comparison of directivity computed for a 
10x10 planar array. The reference directivity is derived 
using (4.24) while the calculated directivity is from (4.10). 



Table 4.1 Planar Array Data Comparison 



ANTENNA 


REFERENCE 


CALCULATED 


DIMENSIONS 


DIRECTIVITY (D 0 ) 


DIRECTIVITY (D 0 ) 


M = N = 10 






d x = d y = 1/2 m 






e o = 30° 


23.67 dB 


24.07 dB 


O 

II 

U1 

0 






Px = Py = -A"/ (2/2) 







Figure 4.2 and Figure 4.3 are comparisons of array factor 
patterns for a 5x5 planar array with patterns from [Ref. 3] 
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and the calculated patterns. Figure 4.2 is a representation 
of the array factor patterns in the x-z and y-z planes, and 
the scale of this figure has been increased in order to 
increase the clarity of the pattern. Figure 4.3 is the array 
pattern in the plane at <p = 45°, where the maximum beam is 
directed at <p Q = 45°. The array patterns produced by this 
application are identical to the patterns given in [Ref. 3]. 



27 




FIGURE 4.2 Array factor pattern of a 
excited, planar array with d x = d y = 
45°, in the x-z and y-z planes 



5x5 element, equally 
A/2 m, 0 O = 30°, 0 O = 
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FIGURE 4.3 Array factor pattern of a 5x5 element, equally 
excited, planar array with d x = d y = X/2 m, 0 Q = 30°, <p 0 = 
45°, in the plane 0 = 45° 
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V. THE CIRCULAR ARRAY ANTENNA 



Circular arrays, or ring arrays, are antennas in which 
the radiating elements are placed on a circle with no elements 
positioned inside the circle. Circular antenna arrays with 
elements positioned inside the circle are planar arrays, and 
these antennas are addressed in Chapter IV. Circular array 
antenna are used for radio direction finding, air and space 
navigation, underground propagation, radar, and sonar systems 
[Ref. 3: p. 274]. As with linear and planar arrays, circular 
arrays with uniformly spaced, equally excited elements are 
narrowband antenna systems due to the dependance of the number 
of elements (N) and the circumference of the circle to the 
wavelength (1) [Ref. 3: p. 277]. 

The following inputs are required to analyze a circular 
array: 

N = number of isotropic radiating elements 
a = radius of the circle 
f = frequency of interest 
6 0 = direction of main lobe at 0 = 0 Q 
0 O = direction of main lobe at 0 = 0 o 
I 0 = antenna feed current 
Z o = characteristic feed impedance 
The first two inputs are data primarily obtained from visual 
data such as photographs, and the ’ indicates input parameters 
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that are either known or estimated. The planar array geometry 
is illustrated in Figure 5.1. 



Z 

A 




FIGURE 5.1 Geometry of N-element circular array of 
isotropic point sources 
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This Mathcad application is written for N-element 
circular arrays with uniformly spaced, equally excited, 



isotropic radiating elements on the x-y plane. The 



polarization of the array is not calculated since isotropic 
point sources are used, and the antenna input resistance (R jn ) 
is assumed to be equivalent to the radiation resistance (R r ) . 
This assumption results in the computation of ideal gain (G) . 
As with the linear and planar arrays, the radiation pattern of 
a specific element type is neglected in the analysis of the 
circular array since in normal usage it will have little 
effect for a circular array with a large number of elements. 
However, as with both the linear and planar arrays, the user 
of this application should be able to use information 
regarding the type of element and the orientation of the 
elements to narrow the range of possible 0 o 's and <p o ' s. 

The normalized array factor (AF) is calculated for far- 
field observations neglecting mutual coupling between adjacent 
elements with equal amplitude excitation and is [Ref. 3: pp. 
274-278; Ref. 4; 350-354]: 




j [jcasinOcos (<(>-4> n ) + ( dimensionless) (5.1) 



where 




n = 1 , 2 , . . . N (radians) 



(5.2) 



a n = -&asin0 o cos ( 4> 0 - $ n ) ( radians ) 



(5.3) 
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Equation (5.1) is the normalized array factor. Equation (5.2) 
is the angular position of the nth element on the x-y plane. 
Equation (5.3) is the phase excitation of the nth element with 
the position of the main beam directed at 0 = 0 o and <p = <p o , 
and k is the wavenumber (2n/\) and a is the radius of the 
array . 

An observation point is considered to be in the far-field 
if all of the following are satisfied [Ref. 4: pp. 24-25): 

r z 1.6 A ( meters ) (5.4) 

r £ 5 (2a) (meters) (5.5) 

r ^ ^ 2 a ^ ( meters ) (5.6) 

A 

The minimum distance to the far-field is found by taking the 
maximum of (5.4), (5.4), and (5.6), where the diameter of the 

circle is taken to be the maximum dimension of the antenna. 

To determine the directivity (D 0 ) for a circular array, 
the approximate equations of radiation intensity (U(0,0)) for 
arrays and radiated power (P rad ) are used. Radiation intensity 
is [Ref. 3: pp. 229-233): 

C7( 0 , 4 > ) = ( AF ) 2 ( W/ solid angle) (5.7) 

The radiated power is [Ref. 3: pp. 28): 

P iad = f 2n r /2 U ( 6 ' 4 >) sin (0) d0d<t> m (5.8) 

Jo Jo 
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In (5.8), P rad is computed by integrating over all angles above 
the x-y plane (0 < 0 < n/2) [Ref. 7: p. 78]. This assumes 
that the array has a pattern maxima above the x-y plane with 
no radiation below the x-y plane. The directivity for the 
circular array is [Ref. 3: pp. 229-233]: 

4 7t U 

D 0 = — ( dimensionless ) (5.9) 

P iad 

In (5.9) , U maK is unity and occurs at 0 = 0 O and <p = <p 0 since 
the array factor is normalized. 

Effective isotropic radiated power (EIRP) is determined 
by [ Ref . 4 : p . 62 ] : 

EIRP = P zad D 0 (W) (5.10) 

Assuming that the power delivered to a circular array of 
N isotropic, uniformly spaced radiating elements is equivalent 
to the radiated power, we can estimate then the radiation 
resistance (R r ) of the array as [Ref. 3: p. 55]: 

= 2 ! (Q) (5.11) 

Since the radiating elements are assumed to be ideal point 
sources, the input resistance of a circular array (R in ) is 
estimated to be equivalent to the radiation resistance: 

R in = R z (Q) (5.12) 
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Since the input resistance is estimated, the gain (G) of 
a circular array cannot be calculated, but an ideal value can 
be derived. As with linear and planar arrays, addressed in 
Chapter III and Chapter IV, respectively, the total antenna 
efficiency (e t ) is [Ref. 3: pp. 44-45]: 

«c - ■ 1-IIT ( dimensionless ) (5.13) 



where 



e cd = 1 ( dimensionless ) 



(5.14) 



e r = 1 — | r | 2 ( dimensionless ) (5.15) 

r = — — ( dimensionless ) (5.16) 

R in + Z o 

In (5.14), e cd is the conduction and dielectric efficiency, 

which is unity for an ideal lossless antenna. In (5.16), Z 0 

is the characteristic impedance of a transmission line, and T 
is the voltage reflection coefficient. Therefore, the total 
antenna efficiency is equal to the mismatch efficiency of 
(5.15) when an ideal lossless antenna is assumed, and the gain 
of a circular array is (Ref. 3: pp. 43-44]: 

G = t c D a ( dimensionless ) (5.17) 

G(dB) = 10 log 10 (e c D o ) { dB) (5.18) 
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The maximum effective aperture (A ero ) of a circular array 
is derived assuming polarization mismatches are negligible 
(this Mathcad application is constructed using isotropic 
radiating elements for which polarization loss factor (PLF) 
and antenna polarization are not calculated) [Ref. 3: p. 63]: 

= G [ X)2 PLF ( m 2 ) (5.19) 

em 4 m 

PLF = 1 ( dimensionless ) (5.20) 



The effective height (h e J of a circular array is 
estimated by [Ref. 6: p. 42]: 



h em = 2 



^r^em 



N 



(m) 



(5.21) 



In (5.21), r) 0 is the characteristic impedance of free space 
( 12 Ott ) . 

Four two-dimensional array pattern displays are produced 
in this Mathcad application: the first is the x-z plane; the 
second is the y-z plane; the third is a plane perpendicular to 
the x-y plane when 0 = n/4 ; the last is a plane perpendicular 
to the x-y plane when 0 = 0 O . The last display is redundant 
when 0 is either 0, n/4, or n/2. 

The directivity of a uniform circular array approaches 
the number of elements (N) as the radius (a) of the array 
becomes very large [Ref. 3: p. 277]. Table 5.1 provides the 
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directivity of a uniform circular array as a function of the 
radius of the array. 



Table 5.1 Circular Array Directivity Data 



RADIUS (m) 


NUMBER OF 
ELEMENTS 


DIRECTIVITY (D 0 ) 


1.59 


10 


13.70 dB 


10 


10 


12.83 dB 


20 


10 


12.89 dB 


50 


10 


11.14 dB 


100 


10 


10.33 dB 



Figures 5.2 and 5.3 provide comparisons of the array 
radiation patterns computed for a 10-element, uniformly 
spaced, equally excited circular array with a radius of 1.59 
meters with those in [Ref. 3], Figure 5.2 is the radiation 
pattern in the x-z plane, and the radiation pattern in the y-z 
plane is shown in Figure 5.3. For these patterns, 0 O = <p Q = 
0°. The array radiation patterns generated by this 
application are identical to the patterns given in [Ref. 3]. 
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FIGURE 5.2 Array factor pattern of a 10-element, uniformly 
spaced, equally excited circular array with a = 1.59 m, 
and 0 O = 0 O = 0 ° , in the x-z plane 
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FIGURE 5.3 Array factor pattern of a 10-element, uniform 
spaced, equally excited circular array with a = 1.59 m, 
and 0 O = <p 0 = 0°, in the y-z plane 
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VI. THE FOLDED DIPOLE ANTENNA 



The folded dipole is a center-fed wire antenna that 
perforins in a manner similar to that of a linear, center-fed 
thin-wire antenna with the added advantage of a more easily 
matched input impedance [Ref. 8: p. 7-36], The folded dipole 
antenna consists of two parallel dipoles connected at the ends 
to form a narrow wire loop. The folded dipole antenna is 
illustrated in Figure 6.1. 




FIGURE 6.1 Geometry of folded dipole 
antenna 
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The folded dipole can be assumed to behave as a thin-wire 
dipole when the following are satisfied: the length of the 
antenna is much greater than the diameter of the feed 
conductor (L > 2a) ; the distance between the conductors is 
greater than the radius of the feed conductor (d > a) ; 
(27rd/ k) 2 < 1; all conductors are assumed to be perfect [Ref. 
9: p. 172]; 

This Mathcad application calculates the parameters of a 
folded dipole by assuming the antenna behaves as a thin-wire 
dipole with an equivalent radius (a e ) . The current maximum is 
normalized to one amp. The equivalent radius is [Ref. 8: p. 
7-38] : 



a e 



exp 



a 2 ln(a) + ±> 2 ln(£>) + 2abln_(d) 

(a + b ) 2 



(meters) (6.1) 



This application provides the characteristics of the antenna 
both in free space and positioned horizontally over the earth. 
Although the folded dipole positioned vertically over the 
earth occurs frequently in VHF/UHF communications, it is not 
addressed in this research. 

The following inputs are required to analysis a folded 
dipole: 

a = radius of feed conductor 
b = radius of second conductor 
d = spacing between conductors 
L = length of antenna 



41 



h = height of antenna above ground plane 
f = frequency of interest 

e cd = conduction/dielectric efficiency of conductors 
e r = relative dielectric constant of ground plane 
a = conductivity of ground plane 

r ££ = distance for far-field parameter calculations 
a w = incoming wave electric field unit vector for antenna 
in free space 

o a = unit polarization vector for antenna in free space* 
ah w = incoming wave electric field unit vector for 
horizontally positioned antenna 
0h p = antenna polarization direction for horizontally 
positioned antenna 

The first five inputs are physical dimensions obtained through 
photographs, the * indicates input parameters that are either 
known or estimated. 



A. FOLDED DIPOLE IN FREE SPACE 

Assuming the folded dipole behaves as a thin-wire dipole 
with a e < d < L, we have the far-field electric field 
intensity of the antenna aligned with the z-axis as [Ref. 3: 
p. 120]: 



£(0> =j n c 



I„e ~ jkZfc 



cosl— cos (0) \ - cos 



\ 2 



2nr 



ft 



sin (0) 



) - -(f) 



( v/m) 



( 6 . 2 ) 



In (6.2), r] 0 is the intrinsic impedance of free space, equal 
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to 1207 r , and k = 2n/X is the wavenumber. The normalized 
electric field is found by dividing (6.2) by the maximum value 
of the electric field which occurs at 0 = n/2 for L less than 
about 1.251. Clearly, the radiation pattern of the folded 
dipole is identical to that of a conventional dipole of length 
L. 

An observation point is considered to be in the far-field 
if all of the following are satisfied [Ref. 4: pp. 24-25]: 

r £ 1.61 (meters) (6.3) 

r ^ 5 L ( meters ) (6.4) 

r ^ - y (meters) (6.5) 

The minimum distance to the far-field is found by comparing 
the values of (6.3), (6.4), and (6.5) and selecting the 

maximum value. The distance between the conductors is not 
taken into account in determining the minimum distance to the 
far-field since it is assumed that the folded dipole has an 
over all effective radius (a e ) . 

Directivity (D 0 ) for the folded dipole is found by 
determining the radiation intensity (U(0)) and radiated power 
(P rad ) . Since the radiation pattern is symmetric in the x-y 
plane, U(0) and P rad are functions of 0 only. Radiation 
intensity is [Ref. 3: p. 28]: 
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T ^ 

U{ 0 ) = — — | £ 7 ( 0 ) | 2 ( W/ solid angle) ( 6 . 6 ) 

2T 1o 

The radiated power is [Ref. 3: pp. 28]: 

P rad = 2nf n U(6) sin(6) d6 (W) (6.7) 

J o 

Directivity for a folded dipole is [Ref. 3: pp. 29]: 

4 7t U 

D 0 = — ( dimensionless ) ( 6 . 8 ) 

P zad 

where U max = U(tt/2) . 

Effective isotropic radiated power (EIRP) is determined 
by [Ref. 4 : p. 62 ] : 

EIRP = P rad D 0 (W) (6.9) 

The folded dipole is unique in its capability to act as 
a step up impedance transformer when the antenna is resonant. 
The antenna effectively operates as an unbalanced transmission 
line, and it can be analyzed by decomposing the current into 
two distinct modes: a transmission line mode [Figure 6.2] and 
an antenna mode [Figure 6.3]. [Ref. 3: p. 341; Ref. 4: pp. 
205-207] 
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FIGURE 6.2 Transmission 



line mode 




The input impedance (Z jn ) at the feed for the folded 
dipole is given in (6.10), and the radiation resistance (R jn ) 
for the antenna is the real component of (Z jn )) [Ref. 3: 
p.342] : 



Z in = 



f Z ’c 

2Z, + Z 



<Q) 



( 6 . 10 ) 
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In (6.10), Z t is the input impedance of the transmission line 
mode and is obtained from the impedance transfer equation 

Z t =jZ 0 tan|^) (G) (6.11) 



where Z o in (6.11) is the characteristic impedance of a two- 
wire transmission line and is given in (6.12) [Ref. 3: pp. 
341-342; Ref. 8: p. 7-39]: 

Z = — cosh _1 ( — — — ] (Q) (6.12) 

n \2jal 5) 

The input impedance of a linear dipole of length L and 
equivalent radius a e is calculated by [Ref. 3: pp. 124, 127, 
294 , and 342] : 

Z d = R d + jX d (G) (6.13) 



*d = 




(G) 



(6.14) 



*d = 




(G) 



(6.15) 



R 



r 





i- In (kL) - C d (kL)] 

-|-sin(icL) [S i (2kL) - 2S i (kL)] 
cos ( kL) (y + In 



kL 

2 



j + C i ( 2 kL ) -2 C i ( kL ) \ 



(6.16) 
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X, = 



n 0 



-\2SAkL) + cos ( kL ) [2SAkL)-SA2kL)] 

7t l 1 1 1 



(6.17) 



_Ho 
4 TC 



2 if ^ 2 

sin ( kL) \2C i ( kL) -C i ( 2kL) — c^f 



In (6.16) and (6.17), C^x) and S^x) are the cosine and sine 
integrals, respectively, [Ref. 3: pp. 124 and 743-746]: 



C d (x) 



= f J, cos(T) dt = Y.ln(x) - V 
J ~ 1 



( -I)”' 1 * 2 ” 
2n(2n ) ! 



(6.18) 



Sl U) ./ 



x sin ( t ) 



* - E 



n = l 



( ) /i-l ^ (2/3-1) 

(273-1) (2/3-1) ! 



(6.19) 



and y = 0.57721 (Euler's constant). 

Mathcad computes summations more quickly than integrals. 
In this application, the summation forms of (6.18) and (6.19) 
are used with the index n ranging from 1 to 50 to ensure 
proper convergence. The maximum value of the argument of 
either C A (x) or Sj(x) is x = 47TL/X, and for L = X the 

succeeding terms in the summation rapidly decrease for n > 20. 

The gain (G) of a folded dipole is obtained as the 
product of the antenna reflection efficiency (e t ) and the 
directivity of (6.8) where the antenna reflection efficiency 
is [Ref. 3: pp. 43-45]: 

e t = t r t cd = e cd ( 1 — | r | 2 ) ( dimensionless ) (6.20) 
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where 



e r = l - I r I 2 ( dimensionless ) (6.21) 

r = — — (dimensionless) (6.22) 

R in + Z o 

In (6.21) e p is the mismatch efficiency and e cd is the 
conduction and dielectric efficiency in (6.22). In (6.22), Z 0 
is the characteristic impedance of the transmission line given 
in (6.12). Therefore, the gain of the antenna is: 

G = t c D 0 ( dimensionless ) (6.23) 

G(dB) = 10 log 10 (t c D o ) ( dB ) (6.24) 

The folded dipole in free space aligned on the z-axis is 
horizontally polarized with respect to the x-y plane. 
Therefore, given the direction of the incoming wave's unit 
polarization vector (a g ) , the polarization loss factor (PLF) 
is [Ref. 3: pp. 48-53]: 

PLF = \o w ‘ o’ a \ 2 ( dimensionless ) (6.25) 

The maximum effective aperture (A^) of a folded dipole 
is [Ref. 3: p.63]: 

A em = - ( - ) 2 PLF (m 2 ) (6.26) 

4 71 
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The effective height (h em ) of the folded dipole is [Ref. 



6: p. 42] : 



K* = 2 






N 



(m) 



(6.27) 



The bandwidth of the folded dipole antenna is dependant 
on the length of the wires and has better bandwidth 
characteristics than a conventional dipole of the same length. 
It can be assumed that the bandwidth of the folded dipole is 
essentially the same as that of a conventional dipole of an 
radius a e (a < a e < d/2) [Ref. 3: p. 346]. The bandwidth is 
increased when the distance between the conductors (d) is 
increased and/or when the conductor radii a and b are 
increased [Ref. 10: p. 190]. 

Tables 6.1 and 6.2 are comparisons of measured and 
calculated input impedance data for the folded dipole as a 
function of its electrical length. The radii of the 
conductors (a, b) are the same (0.0005 m) while the spacing 
between the conductors (d) is held constant at 0.00625 m and 
A. = 1 m [Ref. 8: pp. 7-40 - 7-41]. 
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TABLE 6.1 Folded Dipole Antenna Data Comparison 



ANTENNA 


MEASURED 


CALCULATED 


LENGTH (L) 


Re(Z jn ) n 


Re(Z jn ) a 


X/n 


2400 


6533 ** 


2X/n 


2400 


2592 


3X/n 


liO 


c— 1 
• 

O 


3A/7T 


300 


324 


5X/7T 


1700 


1685 


6X/n 




1.0 


7 X/n 


100 


100 



50 



TABLE 6.2 Folded Dipole Antenna Data Comparison 



ANTENNA 


MEASURED 


CALCULATED 


LENGTH (L) 


Im(Z in ) n 


Im(Z in ) n 


X/n 


5000 


5437 1 


2.X/ 7T 


-650 


-1204 ** 


AX/n 


-90 


-85 


4 X/tt 


1000 


1061 


5X/TT 


700 


793 


6X/7T 


-190 


-174 


7X/n 


600 


633 



To further test the validity of the Mathcad application, 
a numerical antenna analysis program was also used. The 
program used was ELNEC, a powerful, easy-to-use program for 
modeling and analyzing a wide variety of antenna types 
including ground effects and parasitic structures. The 
fundamental computation portion of ELNEC is the same as 
MININEC (Version 3) , that was developed by the Naval Ocean 
Systems Center [Ref. 11: pp. 5-7]. 

Table 6.3 is a comparison of the gain for a folded dipole 
in free space with different dimensions [Ref. 8: p. 3-29]. 
For each case, the radii of the conductors (0.0005 m) and the 
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distance between the conductors (12.5a m) are the same. As 
can be seen, the comparison between the gains obtained from 
various sources is very good. 



TABLE 6.3 Folded Dipole Antenna Data Comparison 



ANTENNA 


MEASURED 


ELNEC 


Mathcad 


DIMENSION 


Gain (G) 


Gain (G) 


Gain (G) 


L = 1/4 m 


1.7# dB 


1.70 dB 


1.64 dB 


d = L/13 m 








L = 1/2 m 


2.14 dB 


2.20 dB 


2.15 dB 


d - L/ 2 5 . 5 m 









The computed parameters in Tables 6.1, 6.2, and 6.3 are 
essentially the same as those measured or computed with ELNEC 
with exception of the ones marked with **. These discrepancies 
result in using the transmission line model and equivalent 
radius to analysis the folded dipole in this application. 

As previously mentioned, the folded dipole is a center- 
fed wire antenna that performs in a manner similar to a 
linear, center-fed, thin-wire antenna. Therefore, the 
radiation pattern of a folded dipole is also similar to that 
of a thin-wire antenna. The folded dipole radiation pattern 
in free space is illustrated in Figure 6.4. 
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FIGURE 6.4 Radiation Pattern folded dipole in free space 
with L = A/2 m, a = b = 0.0005 m, and a e = 0.004 in 



Figure 6.5 is a comparison of the radiation pattern of a 
folded dipole in free space obtained using ELNEC and the 
Mathcad application. The patterns are not identical, but are 
similar and oriented in the same direction as expected. The 
difference can be attributed to the fact that ELNEC takes into 
account mutual coupling between the conductors while the 
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Mathcad application models the folded dipole as a thin-wire, 
center-fed dipole. 




Figure 6.5 Comparison of radiation patterns using ELNEC 
and Mathcad of a Folded dipole in free space with L = A/2 
m, a = b = 0.0005 m, and a e = 0.004 m 
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B. FOLDED DIPOLE POSITIONED HORIZONTALLY OVER THE EARTH 



Through the use of image theory, the far-field electric 
field intensity of a folded dipole oriented horizontally above 
a flat earth can be obtained. The 0 component of the electric 
field of a horizontal dipole (parallel to the y-axis) in free 
space is [Ref. 3: pp. 143-144]: 



From image theory, the 0 component of the electric field 
over ground is [Ref. 4: pp. 229-235]: 



In (6.28) and (6.29), E 0 (0,0) is valid only above the ground 
(0 < 0 < n/2) [Ref. 3: pp. 135-142]. In (6.29), h is the 
distance from the earth to the antenna and the vertical 
reflection coefficient (T v ) of the ground is [Ref. 4: pp. 231- 
233] : 






>/l - sin 2 (0) sin 2 (<|>) 



£ 0 ^( 0 ) =£ 0 ( 0 , 4 )) (e^ cos(e> - r v e-^ cos(6) ) ( V/m) (6.29) 



e / r cos(0) - Jz' t - sin 2 ( 0 ) 



( dimensionless ) (6.30) 



e^cos^) + yfe 7 ^ - sin 2 ( 0 ) 



where 




O 



( dimens i onl ess) 



(6.31) 
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In (6.30) and (6.31), (e' r ) is the relative complex effective 

dielectric constant of the ground, (e r ) is the relative 
dielectric constant of the ground, and (a) is the conductivity 
of the ground. For a perfect ground plane, o -* <x> and r v = 1. 

The $ component of the electric field of a horizontal 
dipole (parallel to the y-axis) in free space is [Ref. 3: pp. 
135-142 ] : 



£4>(0,<J>) = jn 0 J ° e 



-Jkiff 



2 nr 



££ 



cos|-^sin (0) sin(4>)j -cos| 
\[i - sin 2 (0) sin 2 (<(>) 



kL 

2 



( V/m) 



(6.32) 



Now from image theory, the $ component of the electric field 
over earth is [Ref. 4: pp. 229-235]: 

£4)^(0) =£’4>(0,4>)(e^ cos{0) + r A e-^ cos(e) ) (V/m) (6.33) 



In (6.33), (T h ) is the horizontal reflection coefficient of 
the ground and is [Ref. 4: pp. 231-232]: 



cos(0) - Je' r - sin 2 ( 0 ) 
cos(0) + s ft r z - sin 2 ( 0 ) 



(dimensionless) (6.34) 



When o -* <x>, r h = -1. 

The radiation intensity (U h (0,<£)) of a folded dipole over 
a earth must take into account both the 0 and § components of 
the electric field. Therefore, the radiation intensity is 
[Ref. 3: p. 28]: 
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U h {Q) = I 2 + I 2 ] ( W/ solid angle) (6-35) 



The radiated power (Ph rad ) follows (6.7) except that the field 
is present only above the earth. Thus, 

Jt 

Ph rad = 2n f 2 U h (Q) s in (0)d0 (W) (6.36) 

J o 



Directivity (Dh 0 ) and effective isotropic radiated power 
(EIRPh) for a horizontal folded dipole over the earth are 
identical to (6.8) and (6.9), respectively, with the results 
of (6.36) inserted. 

Input resistance (Rh in ) for the folded dipole over the 
earth plane is found by (where I Q in (6.32) is normalized to 
one amp) [Ref. 3: p. 124]: 




where the radiation resistance (Rh r ) is, 

2 ph 

Rh z = 12* (Q) (6.38) 

\lo\ 

The voltage reflection coefficient (TH) , gain (Gh) , 
maximum effective aperture (Ah em ) , and maximum effective height 
(hh em ) for a horizontal folded dipole over the earth are 
determined using the equations for those parameters for a 
folded dipole in free space, as listed above. 
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The Polarization Loss Factor (PLFh) for the antenna over 
the earth can be determined when the direction of the main 
beam of the polarized antenna (0h p ) is estimated. The antenna 
polarization vector (ah a ) is then calculated, and estimating 
the unit polarization vector of the incoming wave (ah„) , the 
PLF is found from (6.25). 

Table 6.4 is a comparison of measured [Ref. 8: p. 3-29], 
ELNEC, and Mathcad software gains for a horizontal folded 
dipole over a reflecting sheet. The radii of the conductors 
(0.0067 m) are the same. As can be seen, the Mathcad software 
gain is similar to the measured and ELNEC gains. 



TABLE 6.4 Folded Dipole Antenna Data Comparison 



ANTENNA 


MEASURED 


ELNEC 


Mathcad 


DIMENSION 


Gain (G) 


Gci in ( G ) 


Gain (G) 


L = X/2 m 








d = L/25.5 m 


7.14 dB 


6.98 dB 


7.19 dB 


h = X/8 m 









No input impedance data has been found to compare a 
folded dipole oriented horizontally over the earth with that 
computed by this Mathcad application. 

Figures 6.6 and 6.7 are the H-plane and E-plane radiation 
patterns of a folded dipole oriented horizontally over the 
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earth, respectively, 
thin-wire center-fed 
below are similar to 



Since the folded dipole is similar to a 
dipole of length L, the illustrations 
those of a dipole of length L. 




FIGURE 6.6 H-plane radiation pattern of a folded dipole 
with L = X/2 m, a = b = 0.0005 m, d = 0.00625 m, h = 0.5 
m, o = 0.01 S/m, e r = 15 
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FIGURE 6.7 E-plane radiation pattern of a folded dipole 
with L = A./2 m, a = b = 0.0005 m, d = 0.00625 m, h = 0.5 
m, a = 0.01 S/m, e r = 15 



Figures 6.8 and 6.9 are comparisons of the H-plane and E- 
plane radiation patterns between those obtained with ELNEC and 
Mathcad software of a folded dipole oriented horizontally over 
the earth, respectively. As expected, the patterns are 
similar. As previosly mentioned, the difference can be 
attributed to the fact that ELNEC takes into account mutual 
coupling between the conductors while the Mathcad application 
is modeled after a thin-wire, center-fed dipole. 
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Figure 6.8 Comparison between ELNEC and Mathcad of H-plane 
radiation patterns of a folded dipole with L = 1/2 m, a = 
b = 0.0005 m, d = 0.00625 m, h = 0.5 m, a = 0.01 S/m, e r = 
15 




Figure 6.9 Comparison between ELNEC and Mathcad of E-plane 
radiation patterns of a folded dipole with L = 1/2 m, a = 
b = 0.0005 m, d = 0.00625 m, h = 0.5 m, o = 0.01 S/m, e r = 
15 
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VII, THE CAGED DIPOLE ANTENNA 



The caged dipole antenna is used in communications for 
its broadband antenna characteristics. Its performance is 
similar to that of a thick cylindrical dipole. Since thick 
cylindrical dipoles are usually too clumsy to use except at 
very short wavelengths, thick dipoles are replaced with a 
caged dipole which consists of a number of thin wire dipoles 
arranged to form a circular cross section. The antenna is 
usually supported by a central conductor and support members 
that provide rigidity but contribute little to the electric 
properties of the caged dipole. The caged dipole (without the 
central conductor and the support members) is illustrated in 
Figure 7.1 [Ref. 12: pp. 172-174]. 





FIGURE 7 . 1 Geometry of caged dipole antenna 
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This Mathcad application models the caged dipole antenna 
as a circular array of center-fed, equally excited, uniformly 
spaced, thin-wire dipoles. This approach is taken to account 
for the number of conductors that are present in the caged 
dipole. As with the arrays previously considered, mutual 
coupling between the dipoles is neglected. In this manner, 
this application simulates the characteristics of a thick 
cylindrical dipole. The thin-wire dipoles are modeled after 
finite length dipoles in free space that have their length 
much greater than their diameter (L > wire diameter) , and the 
dipoles are assumed to be perfect conductors. Since the 
central conductor and the support wires have little effect on 
the electrical properties of a caged dipole, mutual coupling 
from these components are not taken into account in this 
application. As seen in Figure 7.1, the inner ends of a caged 
dipole, the feed points, maybe coned in order to reduce base 
capacity [Ref. 12: pp. 172-174]. This effect is also not 
addressed . 

This application provides the characteristics of the 
antenna in free space as well as and oriented both vertically 
and horizontally over the earth. 



A. CAGED DIPOLE IN FREE SPACE 

The following inputs are required to analyze a caged 
dipole in free space: 
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N = number of conductors 
a = radius of antenna 
L = length of antenna 
f = frequency of interest 
I Q = antenna feed current 
Z G = input impedance 

e cd = conduction/dielectric efficiency of conductors 
r ff = distance for far-field parameter calculations 
a w = incoming wave electric field unit vector for 
antenna 

a a = unit polarization vector for antenna 
The first three inputs are physical dimensions obtained 
primarily from photographs, and the * indicates input 
parameters that are either known or estimated. 

Assuming the caged dipole behaves as a circular array 
with equally excited (I n ) , uniformly spaced elements in the x- 
y plane with a radius of (a) , we have the array factor 
( AF (0,0) ) [Ref. 3: pp. 274-278]: 

AF ( 0 , 4> ) = £ J ne 7 t^s in 0 c ° s «|,-<p n »] ( dimensionless ) (7.1) 

n = l 

where the angular position in the nth element in the x-y plane 
is : 

= 27t (-^)' n = 1, 2, . . . N (radians) (7.2) 
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Equation (7.1) is derived from: 



N 



AF(0,<t>) = £ Jne j[*asinecos( 4 >-» n )*« n } ( dimensi onl ess) (7.3) 



where the phase excitation of the nth element is: 



a. n = 0 = -Jcasin0 o cos ( <J> 0 - $> n ) ( radians ) 



(7.4) 



In (7.3), and the phase excitation of the nth element (a n ) is 
zero as the result of desiring the main beam of the array to 
be directed perpendicular to the length of the antenna (0 O = 
7r/ 2). In (7.1), k is the wavenumber (2ir/k). 

Assuming each conductor of a caged dipole behaves as a 
thin-wire dipole, we have the far-field electric field 
intensity of a single conductor aligned with the z-axis as 
[Ref. 3: p. 120]: 



In (7.5), r? Q = 120?r and is the intrinsic impedance of free 
space . 

As previously discussed in Chapter III, the total 
electric field of an array constructed with identical elements 
is the product of the field of a single element and the array 
factor as a consequence of the principal of pattern 
multiplication [Ref. 3: p. 207]. Therefore, the total 
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electric field of a caged dipole in free space is approximated 
as : 

£ t (0,4>) = E(B) AF(0,4>) (V/m) (7.6) 

An observation point is considered to be in the far-field 
if all of the following are satisfied [Ref. 4: pp. 24-25]: 

r £ 1.6X ( meters ) (7.7) 



r ^ 5 s/T? + (2 a) 2 



r ^ 2 (L 2 + (2a) 2 ) 



(meters) 



( meters ) 



(7.8) 



(7.9) 



The minimum distance to the far-field is found by taking the 
maximum of (7.7) through (7.9). In (7.8) and (7.9), the 
maximum dimension of the antenna is the diagonal length of the 
array. 

Directivity (D 0 ) for a caged dipole is found by 
determining the radiation intensity (U(0,0)) and radiated 
power (P rad ) . Radiation intensity is [Ref. 3: p. 28]: 

r* 2 

(7(0,4)) = — — \E t ( 0, 4> ) | 2 ( w/ solid angle) (7.10) 

2 ti 0 

The radiated power is [Ref. 3: pp. 28]: 

P rad = f o 2n fJu(e,<t>) sin(O) d6d<t> (HO (7.11) 
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Directivity for a caged dipole is [Ref. 3: pp. 29]: 

4 jt U 

D 0 = — (dimensionless) (7.12) 

P iad 

where U max = U(tt/2,0) . 

Effective isotropic radiated power (EIRP) is determined 
by [Ref. 4 : p. 62 ] : 

EIRP = P iad D 0 (W) (7.13) 



Input resistance (R in ) for a caged dipole in free space 
is found from [Ref. 3: p. 124]: 

R r 



*in = 



sin 



If 



(Cl) 



(7.14) 



where the radiation resistance (R r ) is: 

2P r 



R r = 



rad 



1 \I„\ 



(Cl) 



(7.15) 



The gain (G) of a caged dipole is obtained as the product 
of the antenna reflection efficiency (e t ) and the directivity, 
where the antenna reflection efficiency is ([Ref. 3: pp. 43- 
45] : 

e t = t z t cd = e cd (l-|r| 2 ) ( dimensionless ) (7.16) 



where 

e r = 1 - | T | 2 (dimensionless) (7.17) 
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( dimensionless ) 



(7.18) 



p = ^in 

Rin + Z° 



The parameter e r is the mismatch efficiency, e cd is the 
conduction and dielectric efficiency which is unity for a 
lossless conductor, and Z 0 is the characteristic impedance of 
a transmission line. Therefore, the gain of the antenna is: 

G = t c D a ( dimensionless ) (7.19) 

G(dB) = 10 log 10 (e c D o ) ( dB ) (7.20) 

A caged dipole in free space aligned parallel to the z- 
axis is horizontally polarized with respect to the x-y plane. 
Therefore, given the direction of the incoming wave's unit 
polarization vector (a w ) , the polarization loss factor (PLF) 
is [Ref. 3: pp. 48-53]: 

PLF = |S w -o’| 2 ( dimensionless ) (7.21) 

The maximum effective aperture (A em ) of a caged dipole is 
[Ref. 3: p.63]: 

A em = G( / )2 PLF ( m 2 ) (7.22) 

em 4 k 



The effective height (h em ) of the caged dipole is [Ref. 
6 : p . 42 ] : 



= 2 



R i A em 



N 



im) 



(7.23) 
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The bandwidth of a thick cylindrical dipole antenna is 
dependent on the ratio of the length of the antenna to the 
diameter of the antenna (L/2a) . When L/2a ~ 5000 the 

acceptable bandwidth of a thick cylindrical dipole is about 3% 
of the center frequency, and when L/2a « 260 the bandwidth is 
about 30%. Hence, if we assume a linear relationship between 
bandwidth and L/2a, and since the caged dipole's performance 
is similar to that of a thick cylindrical dipole, then the 
bandwidth of a caged dipole is [Ref. 3: p. 333]: 

BN*f high -f lw (HZ) (7.24) 

where the upper and lower frequencies are: 

W = iHZ) <7 - 25) 



f 1 ow ^ c 




(Hz) 



(7.26) 



In (7.25) and (7.25), f c is the center frequency and F(L/2a) 
is the linear relationship function between antenna length and 
diameter. This relationship is obtained by solving y = mx + 
f c for m using Xj = 5000, y l = 3% and x 2 = 260, y 2 = 30%. The 
relationship becomes: 

f( — U -5 .696-10* 5 /— Uo.3148 ( dimensionless ) (7.27) 
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Table 7.1 is a comparison of the gain of a thick 
cylindrical dipole with that of a caged dipole (N = 8) in free 
space computed both by ELNEC and Mathcad software [Ref. 8: p. 
2-38]. As previously mentioned in Chapter VII, ELNEC is an 
easy-to-use numerical analysis program for modeling and 
analyzing different antenna types [Ref. 11: pp. 5-7]. As can 
be seen, the comparison between the measured, ELNEC, and 
Mathcad application values are similar. 



TABLE 7.1 Caged Dipole Antenna Data Comparison 



ANTENNA 


MEASURED 


ELNEC 


Mathcad 


DIMENSIONS 


Gain (G) 


Go. in ( G ) 


Gain (G) 


L = 1/2 m 








a = 0.025 m 


2.14 dB 


1.87 dB 


1.88 dB 


N = 8 

1- - 









Figure 7.2 is the radiation pattern of a caged dipole (N 
= 8) in free space. It is essentially equivalent to that of 
a center-fed, thick cylindrical dipole having the same 
conductor length and diameter [Ref. 8: p. 2-38]: 
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FIGURE 7.2 Radiation pattern of a caged dipole in free 
space with L = X/2 m, N = 8, and a = 0.125 m 



Figure 7.3 is a comparison of the radiation pattern of a 
caged dipole (N = 8) obtained using ELNEC and Mathcad. The 
patterns are not identical, but are similar and oriented in 
the same direction as expected. The difference can be 
attributed to the fact that ELNEC takes into account mutual 
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coupling between elements while the Mathcad application 



ignores these effects. 




FIGURE 7.3 Comparison of radiation patterns using ELNEC 
and Mathcad of a caged dipole in free space with L = A./2 
m, N = 8, and a = 0.125 m 
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B. CAGED DIPOLE ORIENTED VERTICALLY OVER EARTH 



The following known or estimated inputs are additional 
data required to analyze a caged dipole oriented vertically 
over the earth. The height of the antenna above ground is a 
physical dimension obtained through photographs: 
h = antenna height above ground 
e r = relative dielectric constant of earth 
o = conductivity of earth 

av w = incoming wave electric field unit vector for 
vertical antenna 

0v p = antenna polarization direction for vertically 
orientated antenna 

As with the caged dipole in free space, the analysis of 
a caged dipole positioned vertically above the ground does not 
take into account the central conductor and support members 
since they produce negligible effects on the antenna's 
electrical properties. 

From image theory and pattern multiplication for arrays, 
the total far-field electric field intensity of a caged dipole 
oriented vertically above a flat earth is [Ref. 4: pp. 229- 
235] : 

£V t (e,4>) =Ev(Q) AFv(Q,4>){e jkhcos{e) + r v e-i khcosW ) (V/m) (7.28) 

In (7.28), Ev(0) is the far-field electric field intensity of 
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a finite length dipole in free space which is given by (7.5), 
and AFv(0,0) is the array factor for a circular array which is 
(7.1). The total electric field intensity is valid only above 
the ground (0 < 0 < n/2) [Ref. 3: pp. 135-142]. The distance 
from the earth to the center of the antenna is (h) , and the 
vertical reflection coefficient ( F v ) of the ground is [Ref. 4: 
pp. 231- 233 ] : 



r v (0) 



e^cos^) - Je7~ z - sin 2 (0) 
e^cosf©) + ^/e^-sin 2 (0) 



( dimensionless ) (7.29) 



where 

z‘ r = e - j — ( dimensionless ) (7.30) 

r r 2 tc f e 0 v ' 

In (7.29) and (7.30), (e' r ) is the relative complex effective 
dielectric constant of the ground, (e r ) is the relative 
dielectric constant of the ground, (a) is the conductivity of 
the ground, and e G is the permittivity of free space (l/(367r) 
x 10~ 9 ) . For a perfect ground plane, a -*■ « and T v = 1 . 

The radiation intensity (U v (0,0)) of a caged dipole over 
earth is [Ref. 3: p. 28]: 

r* - 

t/ v (0, <[>) = — —[\Ev t (0, (j>) | 2 ] (W/ solid angle) (7.31) 



The radiated power (Pv rad ) is present only above the earth. 
Thus, [Ref. 3 : p. 28] : 



(7.32) 



Pv rad = f 2 7 2 ^(e / <D)sin(e ) d0d<D (W) 

Directivity (Dv 0 ) and effective isotropic radiated power 
(EIRPv) for a vertically oriented caged dipole over the earth 
are given by (7.12) and (7.13), respectively, with the results 
of (7.32) inserted. Uv max in (7.12) is simply the maximum 
value of (7.31) which, in Mathcad, is found by constructing 
(7.31) as a vector and then finding the maximum value of the 
vector . 

Input resistance (Rv in ) , voltage reflection coefficient 
(TV) , gain (Gh) , maximum effective aperture (Av em ) , and maximum 
effective height (hv em ) for a vertical caged dipole over the 
earth are determined using the equations for those parameters 
for a caged dipole in free space. 

The Polarization Loss Factor (PLFv) for the antenna over 
the earth can be determined when the direction of the main 
beam of the polarized antenna (0v p ) is estimated. The antenna 
polarization vector (av a ) is then calculated and, estimating 
the unit polarization vector of the incoming wave (crv w ) , PLF 
is found from (7.21). 

Table 7.2 is a comparison of the gain of a vertical caged 
dipole over a flat earth obtained with both ELNEC and the 
Mathcad application. As can be seen, the gains are similar. 
As previously mentioned, the difference in the gains can be 
attributed to the fact that ELNEC takes into account mutual 
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coupling between the elements as well as ground effects. The 
Mathcad application ignores mutual coupling between elements 
and approximates ground effects. 



TABLE 7.2 Caged Dipole Antenna Gain Comparison 



ANTENNA 

DIMENSIONS 


ELNEC 
Gci in ( G ) 


Mathcad 
Gain (G) 


L = 1/2 m 
a = 0.125 m 






N = 8 


4.25 dB 


4.50 dB 


h = 0.25 m 
a = 0.01 S/m 






e r = 15 







Since the caged dipole behaves in a manner similar to 
that of a thin dipole except for bandwidth characteristics, 
the radiation pattern of a vertically oriented caged dipole 
over ground will resemble that of a thin vertical dipole over 
ground. Figure 7.4 is a plot of the radiation pattern of a 
caged dipole over earth. There is a very close resemblance to 
the radiation pattern of a thin vertical dipole over ground 
[Ref. 4: p. 234]. Figure 7.5 is a comparison of radiation 
patterns computed by ELNEC and Mathcad software. As expected, 
the patterns are similar. 
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FIGURE 7.4 E-plane radiation pattern of a caged dipole 
with L = X/2 m, N = 8, a = 0.125 m, h = 0.25 m, a - 0.01 

S/m, and e r = 15 




— Wathcad Pattern. ... ELNEC Pattern 



FIGURE 7.5 E-plane radiation pattern comparison between 
ELNEC and Mathcad of a caged dipole with L = X/2 m, N = 
8, a = 0.125 m, h = 0.25 m, a = 0.01 S/m, and e r = 15 
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C. CAGED DIPOLE ORIENTED HORIZONTALLY OVER EARTH 



The following known or estimated inputs are additional 
data required to analyze a caged dipole oriented horizontally 
over the earth: 

h = antenna height above ground 

e r = relative dielectric constant of earth 

a = conductivity of earth 

crh w = incoming wave electric field unit vector for 
vertical antenna 

0h p = antenna polarization direction for vertically 
orientated antenna 

As with the evaluations of a caged dipole both in free 
space and oriented vertically over the earth, the analysis of 
a caged dipole positioned horizontally above the ground does 
not take into account the central conductor and support 
members since they produce negligible effects on the antenna's 
electrical properties. 

Through the use of image theory, the far-field electric 
field intensity of a caged dipole oriented horizontally above 
a flat earth can be obtained. The 0 component of the electric 
field of a horizontal dipole (parallel to the y-axis) in free 
space is [Ref. 3 : pp. 143 - 144 ]: 
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£8A(9,<t>) “ill, 



I.e- 3k, « 



2 mr 



ff 



cos|-^sin (0) sin((|>)J - cos|-^ 
y/l - sin 2 (0) sin 2 (<t>) 



(V/m) 



(7.33) 



The 0 component of the total electric field over ground 
is [Ref. 4: pp. 229-235]: 

EQh c (0, 4>) =£0h(0 / <|>) AFh(Q,<b)(e jkhcos(d) - r v e- jkhcosW ) (v/m) (7.34) 



The array factor (AFh(0,0)) in (7.34) for a horizontal array 
is determined by [Ref. 3: pp. 274 and 776]: 

AFh(Q,<t>) = £ J n e jJcasinesin *sino n ( dimensionless ) (7.35) 

n = l 



In (7.33) through (7.35), the electric field (E0h(0,0)), the 
total electric field (E0h t (0 , <p) ) , and array factor (AFh(0,0)) 
are valid only above the ground (0 < 0 < n/2) [Ref. 3: pp. 
135-142]. In (7.34), (h) is the distance from the earth to 
the center of the antenna, and the vertical reflection 
coefficient (T v ) of the ground is given by (7.29). 

The $ component of the electric field of a horizontal 
dipole (parallel to the y-axis) in free space is [Ref. 3: pp. 
135-142] : 



I e ~ JkZ(f 

Etyh (0, <J>) =jr\ 0 n 

2TZZ ff 



cosj^sin (0) sin(<(>)J 



—Hr) 



y/l - sin 2 (0) sin 2 (4>) 



( V/m) 



(7.36) 



From image theory, the $ component of the total electric field 
over earth is [Ref. 4: pp. 229-235]: 
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E$h t (0, <J>) =E<\>h (0, <t>) AFh (0,0) (e jA: * cos(0) + r h e-^ hcos(0) ) (V/m) (7.37) 



In (7.37), (AFh (0 , 0) ) is given by (7.35), and (T h ) is the 
horizontal reflection coefficient of the ground and is [Ref. 
4: pp. 231-232]: 



cos ( 0 ) - Jz' T - sin 2 ( 0 ) 
cos(0) + r - sin 2 ( 0 ) 



(dimensionless) (7.38) 



When a - oo, r h = -1. 

The radiation intensity (U h (0,0)) of a caged dipole over 
a earth must take into account both the 0 and $ components of 
the electric field. Therefore, the radiation intensity is 
[Ref. 3: p. 28]: 

r* 2 

Uh(Q,<\>) =_IU|£eh t (0,<i>) | 2+ |£0A t (0,0) | 2 1 (W/ solid angle) (7.39) 

2t| 0 1 j 

The radiated power (Ph rad ) is present only above the earth. 
Thus, 

ph rad = [ Z * f" /2 Uh(Q,<t>) sin(0) d0d<}) ( W) (7.40) 

Jo Jo 



The remaining parameters (i.e., Dh 0 , EIRPh, Rh in , etc.) 
for a horizontal caged dipole over the earth are calculated 
using the equations for these parameters listed previously. 

Table 7.3 is a comparison of the gain of a horizontal 
caged dipole over a flat earth obtained with ELNEC and the 
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Mathcad application. As can be seen, the gains are similar. 
As previously mentioned, the difference in gains can be 
attributed to the fact that ELNEC takes into account mutual 
coupling between the elements as well as ground effects. The 
Mathcad application ignores mutual coupling between elements 
and approximates ground effects. 



TABLE 7.3 Caged Dipole Antenna Gain Comparison 



ANTENNA 

DIMENSIONS 


ELNEC 
Gci in ( G ) 


Mathcad 
Gain (G) 


L = A/2 m 
a = 0.025 m 






N = 8 
h = 1.0 m 
o = 0.01 S/m 
e r = 15 

1 


4.99 dB 


4.73 dB 



As previously mentioned with regard to the vertically 
oriented caged dipole over the earth, the caged dipole behaves 
in a manner similar to that of a thin cylindrical dipole. 
Therefore, the radiation patterns of a horizontally oriented 
caged dipole are similar to those of horizontal thin 
cylindrical dipole which in turn resembles that of a thin 
dipole. Figures 7.6 and 7.7 illustrate the H-plane and E- 
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plane patterns, respectively, of a caged dipole over the earth 
[Ref. 4: pp. 231-237]. The radiation patterns of the 
horizontally orientated caged dipole over a flat earth are 
very similar to those of a horizontal thin dipole over earth. 
Figures 7.8 and 7.9 compare the radiation patterns between 
those computed with ELNEC and those computed by Mathcad 
software. As expected, the patterns are similar. 




FIGURE 7.6 H-plane radiation pattern of a caged dipole 
with L = X/2 m, N = 8, a = 0.125 m, h = 1.0 m, a - 0.01 
S/m, and e r = 15 
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FIGURE 7.7 E-plane radiation pattern of a caged dipole 
with L = A/2 m, N = 8, a = 0.125 m, h = 1.0 m, a = 0.01 
S/m, and e r = 15 




— Mathcod Pattern. ... ELNEC Pattern 



FIGURE 7.8 Comparison between ELNEC and Mathcad of H-plane 
radiation patterns of a caged dipole with L = A/2 m, N = 
8, a = 0.125 m, h = 1.0 m, a = 0.01 S/m, and e r - 15 
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1 




— Mathcad Pattern, ... ELNEC Pattern 



FIGURE 7.9 Comparison between ELNEC and Mathcad of E-plane 
radiation patterns of a caged dipole with L = A/2 m, N = 
8, a = 0.125 m, h = 1.0 m, a = 0.01 S/m, and e r = 15 
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VIII. PARABOLIC REFLECTORS 



Many communications and radar operations require very 
large values of antenna gain (G) that are difficult to achieve 
with a single device. One of the most popular ways to 
construct a very high gain antenna is to reflect the 
electromagnetic energy of a relatively small antenna off a 
metallic parabolic dish in the far-field of the feed device 
[Ref. 13: p. 12-2]. If the source antenna is located at the 
focus of the parabola and if the radius (a) of the parabolic 
dish periphery is large with respect to wavelength (X) , then 
the reflector may be approximated as a uniform phase aperture 
antenna. Gain in excess of 30 dB is common for this type of 
structure [Ref. 4: pp. 423-424]. The geometry of a parabolic 
reflector is illustrated in Figure 8.1, where the parameters 
include the focal distance (f), the distance (p) from the 
origin to any point (P) on the reflector surface, and the 
distance (r') from the origin to the projection of point (P) 
onto the z = 0 plane. The primed angles in Figure 8.1 
correspond to the parameters associated with the feed horn. 
The Mathcad application user should note that the focal 
distance should exceed the minimum distance to far-field for 
the feed antenna. 



85 




The Mathcad parabolic reflector applications are unique 
in that the software is not written as a stand alone package. 
Rather, the reflector software is run as an addendum to the 
programs that already exist for various antennas [Ref. 2: pp. 
4-15, 67-114]. Four types of feeds are analyzed: helical, 
planar equiangular spiral, pyramidal horn, and conical horn. 
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The following additional inputs are required to analyze 
parabolic reflectors: 

foe = reflector focal length 

a ref = radius of the mouth of the reflector 

E 0 = electric field scale factor 

e b = blockage efficiency 

e sp = spar efficiency 

r ££2 = secondary far field observation distance 
x,y,z = polarization loss factor coordinates 
a w = incoming wave unit polarization vector 
tl, t2 = radiated power increments for <p and 0 
il = secondary field increments 
N,M = summation increments 

The first two inputs are physical dimensions obtained 
primarily through photographs, inputs with the * are 
parameters that are either known or estimated, and the last 
inputs are only used to affect the processing time of the 
application. 

In order to calculate the far-field patterns of a 
parabolic reflector, the electric field at the aperture of the 
device must be known. Since the electric field patterns of 
interest are in the far-field, it is assumed that the circular 
aperture created by the reflector is located at the origin and 
is parallel to the z = 0 plane. The Cartesian components of 
the aperture electric field are expressed as [Ref. 14: pp. 

121-123; Ref. 15: pp. 29-31]: 
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E a*( r ' >$') = [- E w (6' (r') ,0 / )cos(0 / ) 

+ E, f (0' (r 7 ) ,<p') sin(0 7 ) ]e 



jk ( p( r' ) -2foc ) 



(V/m) 



( 8 . 1 ) 



E ay (r' ,<p' ) = [E M (0'(r') ,<p')sin(<p') 

+ E^ f (9 / (r / ) , <p' ) cos (<p f ) ]e 



jk ( />( r' ) -2foc ) 



(V/m) 



( 8 . 2 ) 



In (8.1) and (8.2), (E ax (r' ,6' )) and (E ay (r' ,0')) are the 
aperture electric fields and k is the wavenumber ( 2n/X ). 

The relationships between (0' ), (p), and (r' ) needed to 

evaluate (8.1) and (8.2) are expressed as [Ref. 4: p. 426]: 



The electric field from the entire parabolic assembly is 
commonly referred to as the secondary pattern. Referring to 
Figure 8.1, the electric field in the -z half-plane is 
considered the secondary pattern. The secondary field of a 
parabolic reflector is calculated using the electric field 
integral equation solution and computing the component vector 
potentials (P x (0,0)), (P y (0,0)) as follows [Ref. 4: p. 426]: 



0 ; = 2 tan' 1 — ^ — ( radians ) 

^ 2 foe) 



(8.3) 




(m) 



(8.4) 




’ (r / ((/) ^ jfto’ / sin(0)(cos (<t>') cos (4>) ♦ sin(4/) sin(<|)) 



(8.5) 



* r' dr' d 4/ 



(V/m) 
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Py(6,0) = r 2 ” r aref E ay (r' ,0') [ e j kr, sin<8)(' 



cos (<£' ) cos (0 ) ♦ sin(0' ) sin(0) 



>] 



( 8 . 6 ) 



* r' dr ' dtp' 



(V/m) 



Double integrals in Mathcad take significantly longer to 
converge to an answer than do single integrals [Ref. 15: p. 
200]. To reduce the processing time of the applications, 
(8.5) and (8.6) were each separated into two single integrals. 
Therefore, for the applications, (8.5) and (8.6) are each 
computed as: 



AO,^,^ 7 ) = f E a (r', 4 /) e ^ r ' sill(e, [ cos(| l ,/|co s (| l , > + si n «|. / )s in( 4>)] r / c?r / (8.8) 



Any contribution from feed system back lobes is ignored, 
and the component vector potentials calculated in (8.7) and 
(8.8) are applied to the subsequent formulas to compute the 
secondary field [Ref. 4: p. 383]: 




(8.7) 



where 




[P x (0, 4>) cos (<t>) + P y (0, <|> ) sin (<j>) ] (V/m) (8.9) 
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The radiation intensity (U(0,0)), radiated power (P rad ), 
and directivity (D 0 ) of the secondary field pattern are 
expressed using the intrinsic impedance of free space (r) 0 = 
1207 r) as follows [Ref. 3: pp. 27-3.0]: 

U(Q, 4>) = -^^—[\E e (0, 4>) | 2 + (0, 4>) I 2 ] (W/ solid angle) (8.11) 

p rad = f 2 n f n/ 2 c/(e,<t>)sin(e)ddd<t> (w) (8.12) 

Jo Jo 



4ti£7 

^ max 

0 P 

^ 1 ad 



( dimensionless) 



(8.13) 



In (8.13), ( U raax ) is the maximum value of radiation intensity, 
obtained by creating (U max ) as a matrix and finding the maximum 
value of the matrix. 

The computer processing time required to evaluate (8.12) 
is prohibitively long. Thus, the parabolic reflector Mathcad 
applications use the following numerical approximation to 
evaluate radiated power [Ref. 3: pp. 37-42]: 



rad 






(W) 



(8.14) 



In (8.14), (N and M) are user selected integers. This 
approximation for P rad reduces the processing time of the 
applications substantially. As N and M increase, the results 
obtained with (8.14) converge to that obtained with (8.12). 



90 



When N = M = 3, the processing time for only P rad of the 
helical feed reflector option is approximately four hours on 
a 33 MHz, 386 machine, and when N and M is increased to 5, the 
time to calculate P rad for the same application is ten hours. 
For the applications, we chose N = M = 10 to evaluate P rad with 
reasonable accuracy. With this choice, it takes approximately 
one day for the application to run, but this is much less than 
the thirty-six hours it takes if (8.12) is used. 

Parabolic reflectors are members of the aperture antenna 
family. The gain of an aperture antenna is the product of its 
aperture efficiency (e ap ) and directivity. Therefore, the gain 
of a parabolic reflector is written as [Ref. 6: p. 575]: 

G = t ap D 0 ( dimensionless ) (8.15) 

Aperture efficiency is the product of several separate 
terms, including spillover efficiency (c s ) , taper efficiency 
(e t ) , phase (or random surface error) efficiency (e p ) , 
polarization efficiency (e x ) , blockage efficiency (e b ) , spar 
efficiency (e sp ) , and ohmic efficiency (e ohmlc ) . 

Spillover efficiency (e s ) is the fraction of power 
radiated by the feed system that is intercepted and reflected 
by the parabolic dish. Spillover efficiency is the most 
important efficiency term for any feed system and is defined 
as [Ref. 6: pp. 583]: 
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( dimensionless ) (8.16) 



r 2n [°°{j iq/) S in (60 dO'cty' 
_ J 0 Jo 

f 2 ” ["u f (00 sin (00 cf0 x c?4> / 
Jo Jo 



In (8.16), (U f (0 7 )) is the radiation intensity function of the 
feed antenna and (0 O ) is the value of (8.3) when (r 7 ) equals 
the radius of the mouth of the reflector (a re£ ) . 

Taper efficiency (e t ) accounts for the lack of amplitude 
uniformity of the feed pattern on the surface of the 
reflector. Taper efficiency, which is a dimensionless 
quantity, is expressed as [Ref. 3: pp. 626-627]: 



e c = 2cot‘ 



0, 



/>2n /• 0 
JO J C 



0 \ 



4 rc £7,(00 



• tan 



lad 



dB'dV 



(8.17) 



2 n re 0 4nU f (0 ; 



r 7 

JO J C 



sin (00 dd'dfy* 



rad 



Phase efficiency (e p ) is a function of the power loss 
that occurs if the field at the mouth of the dish is not in 
phase at every point in the aperture. Phase efficiency is 
approximately given by [Ref. 4: pp. 434-435]: 



-4H6-10 

e p » e A ( dimensionless ) 



(8.18) 



The power that is lost due to cross-polarized electric 
fields in the aperture of the reflector determines 
polarization efficiency (e x ) . Polarization efficiency is 
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difficult to precisely calculate, but can be reasonably 
approximated as 0.98 [Ref, 4: p. 435]. 

Blockage (e b ) and spar (e sp ) efficiencies are determined 
by the presence of feed structures, support struts, and 
possible signal processing equipment in front of the 
reflector's aperture. Precise computation of blockage and 
spar efficiencies requires an extensive method of moments 
simulation and is not attempted by the Mathcad parabolic 
reflector applications. However, if the radius (a f ) of the 
blocking structure at the focal point and the number of spars 
(N) that are (A/2) thick is known, blockage efficiency and 
spar efficiency can be estimated using Tables 10.1 and 10.2, 
respectively [Ref. 4: p.436]. 



TABLE 8.1 Blockage Efficiency 



3f/c*ref 


0.05 


0.10 


0.20 




0.99 


0.96 


0.84 



TABLE 8.2 Spar Efficiency 





^ref 5/A. 


a ref = 50/A 


a ref = 100 /A 


N = 3 


0.95 


0.99 


1.00 


N = 4 


0 . 94 


0.99 


1.00 
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Ohmic efficiency (e ohraic ) for a reflector is usually very 
high but difficult to accurately calculate. The Mathcad 
parabolic reflector applications assume ohmic efficiency is 
0.98. 

With all reflector efficiency terms known or estimated, 
the aperture efficiency (e ap ) , effective aperture (A e ) , and 
isotropic radiated power (EIRP) are given by [Ref. 3: pp. 623- 
630] : 



^ ap ^ C ^ sp^ ohmic 


( dimens i onl ess) 


(8.19) 




(m 2 ) 


(8.20) 


EIRP = P rad D 0 


(W) 


(8.21) 



The magnitude of the feed current (I Q ) for the source 
antenna cannot be estimated based on measured geometry. 
However, if the feed current is normalized to 1 amp, radiation 
resistance (R r ) for the entire reflector structure is: 

2 P 

Rr = -r-rfj (0) (8.22) 

l-U 

In (8.22), (P rad ) is the radiated power from the aperture of 
the reflector, not of the feed antenna. 

The Mathcad parabolic reflector applications use the 
radiation resistance computed for the parabolic reflector 
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antenna system to determine effective height (h era ) by [Ref. 6: 
p. 42]: 




N 



(8.23) 



The bandwidth of a parabolic reflector is assumed to be 
the same as its feed antenna. 

The polarization loss factor (PLF) of the parabolic 
reflector for a user selected point in the secondary far-field 
and a specific wave unit polarization vector (u w ) is found 
using the electric field components of (8.9) and (8.10) as 
follows [Ref. 3: 51-53]: 



and the Cartesian components of the electric field are [Ref. 
3: p. 776]: 



PLF = \d w 'd* a \ 2 ( dimensionless ) 



(8.24) 



where the unit polarization vector of the antenna is: 



5 a (x,y, z) = 



3. X E X + 3. y Ey+ 3. Z E Z 



X y y 



( dimensionless ) (8.25) 



J\E(x,y, z) | 



E x = E q cos (0) cos (4>) - E (J) sin(4)) ( V/m ) 



(8.26) 



E y = .Egcos (0) sin (<t>) + E^cos ((J>) ( V/m ) 



(8.27) 



E z = -E q sin(0) (V/m) 



(8.28) 
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All of the secondary field parameters calculated by the 
Mathcad parabolic reflector applications assume that the 
observation point is in the far-field. Thus, all of the 
subsequent formulas must hold for the application calculations 
to be valid [Ref. 4: pp. 24-25]: 

r £ 1 . 6 A. (m) (8.29) 

r ^ 10 a zef (m) (8.30) 

r ^ (dimensionless) (8.31) 

A 

Table 8.3 is a comparison of reflector directivity 
computed by the Mathcad applications to predictions made for 
an identical antenna with unity aperture efficiency [Ref. 3: 
p. 634]. As can be seen, the comparisons vary by 

approximately 8 dB. This can be attributed to the fact that 
P rad was approximated using (8.14) instead of (8.12), and also 
to the fact that the numerical tolerance variable (TOL in 
Mathcad's math menu) for the application was increased to 0.5 
verses 0.001 in order to reduce processing time. This 
increase in TOL reduces both the accuracy of (8.14) and the 
processing time of the application. 
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TABLE 8.3 Parabolic Reflector Data Comparison 



ANTENNA 


PREDICTED 


CALCULATED 


PARAMETER 


DATA 


DATA 


DIRECTIVITY 


42 dB 


50.2 dB 



Figure 8.2 is a comparison of measured data fed by an 
optimum horn with that calculated by the Mathcad application 
fed by an equiangular spiral for a parabolic reflector (a ref = 
0.320 m, f = 0.147 m) operating at 35 GHz [Ref. 14: p. 126]. 
As can be seen, the comparison is good. The Mathcad pattern 
of Figure 8.2 is not symmetric which can be attributed to the 
number of increments used to increase the speed of processing. 

The reflectors fed by pyramidal and conical horn antennas 
are not included in this report. The applications have been 
developed and are included on the diskette containing the 
software. The complexity of the software results in 
prohibitively long processing time for the two applications 
(i.e., the pyramidal horn fed reflector application was still 
calculating P rad2 after seven days when run on a 50 MHz, 486 
personal computer using Mathcad 4.0). 



97 




FIGURE 8.2 Comparison of Parabolic Reflector E-Plane 
Electric Field Patterns 
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IX. REMARKS AND CONCLUSIONS 



This thesis and the accompanying software applications 
are intended to provide the Naval Maritime Intelligence Center 
(NAVMARINTCEN) with software for the analysis of various 
antenna types where information regarding the antenna is 
restricted primarily to the antenna type and the physical 
dimensions of the antenna. The goal of this research was to 
provide user friendly software that an engineer familiar with 
basic antenna types can easily use and interpret. Although 
the software was developed for that purpose, the complexity of 
formulas that apply to many of the antenna types of necessity 
reduced the simplicity of the applications. 

The software applications developed are written in 
Mathcad 3 . 1 and are compatible with any IBM personal computer 
that supports Mathcad 3 . 1 or Mathcad 4.0 for Windows. Many of 
the applications require extensive processing time. In order 
to reduce processing time, which results from the complexity 
of various integrals, approximate numerical techniques were 
implemented to evaluate the integrals and the numerical 
tolerance of the more complicated applications was reduced. 
Even with the processing time reduced, some programs require 
up to five days to complete when run with Mathcad 4.0 on a 50 
MHz, 486 personal computer. Applications that did not exceed 
a day to complete were run on a 33 MHz, 386 personal computer 
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with Mathcad 4.0. All applications were run using Mathcad 4.0 
since Mathcad 4.0 is approximately twice as fast as Mathcad 
3.1. As previously mentioned, the applications were written 
in Mathcad 3.1 due to sponsor requirements. Mathcad 3.1 is 
compatible with Mathcad 4.0, but Mathcad 4 . 0 is not compatible 
with Mathcad 3.1. 

For this thesis, a number of the antennas could not be 
adequately researched through the use of texts. In some 
cases, such as the caged dipole, very little information on 
the antenna could be found in either texts or in professional 
journals. For the caged dipole, the antenna is modeled as a 
circular array of center-fed, equally excited, uniformly 
spaced, thin-wire dipoles to take into account the number of 
conductors that are present in the caged dipole. This model 
provides good results when compared to predicted (the caged 
dipole's performance is similar to that of a thick cylindrical 
dipole) and ELNEC results. 

It was found that Mathcad is not the ideal tool to 
analyze antennas due to its long processing time. An 
additional drawback of Mathcad is that when a function is in 
the beginning of an application and the same function is 
placed later in the application, Mathcad does not store the 
values of the function, but re-evaluates the function. This 
re-evaluation increases processing time significantly, 
especially when the function itself takes a very long time 
(hours vice minutes) to reach a numerical solution. 
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Nevertheless, the Mathcad applications developed in this 
research allow the user to analyze several antenna types 
(i.e., arrays with isotropic point sources, parabolic 
reflectors with the various feed antennas) that could not 
easily be evaluated using some antenna numerical analysis 
programs such as NEC, ELNEC, or MININEC. For the other 
antenna types researched, the antenna numerical analysis 
programs mentioned would be easier to setup, run, and provide 
parameters closer to those of a real antenna. The 
disadvantage of this approach is that it requires an ability 
to evaluate antennas using a program such as ELNEC, while the 
Mathcad applications do not require any advanced programming 
knowledge on the part of the user. 
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LINEAR ARRAY 

MATH CAD SOFTWARE ARRAY_LN . MCD 

This linear array application is based on N isotropic point sources 
uniformly spaced and aligned on the z-axis. Mutual interference between 
adjacent sources are ignored when calculating antenna parameters. Since 
isotropic point sources are used as the radiating elements, polarization 
of the antenna can not be determined. The gain can only be idealized 
when it is assumed that the antenna input resistance is equivalent to 
the radiation resistance. Bandwith calculations are not carried out as 
a result of the narrowband characteristics of arrays, and in large part 
to the dependance inter-element spacing (d) has with the wavelength (1) . 
The following limits should be adhered to when analyzing Broadside, 
End-Fire, and Phased (scanning) antennas: 

0 < 0 O < 7 i (Scanning Antennas) 

0 o = n/2 (Broadside Antennas) 

0 o = 0 or 7i (Ordinary End-Fire Antennas) 

(Note: Mathcad equations can not use symbolic subscripts. Therefore, 

parameters do not contain subscripts (i.e., AF n will be written AFn) . ) 



parameters (Items with * indicate parameters 


that 


can be computed but 


their 


results are idealized as the result of 


using 


isotropic point 


sources as the radiating elements): 






X = 


Wavelength 






r min 


= Minimum Distance to the Far-Field 






P = 


Phase Difference Between Adjacent Elements 




V = 


Array Factor Phase Shift 






AF n 


= Normalized Array Factor 






U(0) 


= Radiation Intensity 






p rad 


= Radiated Power 






Do = 


Directivity 






EIRP 


= Effective Isotropic Radiated Power 






R r = 


Radiation Resistance 






R in 


= Antenna Input Resistance * 






r = 


Voltage Reflection Coefficient * 






e t = 


Antenna Efficiency * 






G = 


Gain * 






n 

rt em 


= Maximum Effective Aperture * 






^em 


= Maximum Effective Height * 
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The following known or estimated data must be entered: 



N = Number of Isotropic Radiating Elements 

f = Frequency of Interest 

d = Distance Between Adjacent Elements 

0 O = Direction of Main Lobe 

I Q = Antenna Feed Current 

Z Q = Characteristic Feed Impedence 



Enter input data here 



N := 10 



(elements ) 



d := 0.25 



f := 3 10 8 (Hz) 



0o := 0 



lo-l.O (A) 



Zo:=75.0 



Define constants and calculate wavelength: 



c:= 2.9979- 10 8 (meters/sec) 



T|0 •= 120-71 



X :=- 



(meters/cycle ) 



k:=L-* 



X= 0.999 



(meters/cycle) 



ccd •= 1 



(meters ) 



( radians ) 



(D) 



(Q) 



(m 1 ) 



(dimensionless ) 
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Calculate linear array parameters: 



Define angular offset 9: 

0 := (o- 10' 4 ) ,(— ) .. (2 71 - 10' 4 ) 



Minimum Distance to Far-Field r m i n : 



r 1 .6* X (meters) 



r x := 5 Nd 
2 (N d) 2 



(meters ) 



(meters ) 



rmin := max(r) (meters ) 



rmin = 12.509 



(meters ) 



Array Factor AF n : 
P := - k d cos(0o) 



y(0) •= k d cos(0) -f P 



AFn(0) 



. i 


1 N A \ 


sin 


!— v ( S)) 


sin 


f— ^(0)1 




u / 



( radians ) 



( radians ) 



( radians ) 



( dimensionless ) 
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Radiation Intensity U: 



U(0) • = (AFn(0)) 2 



(W/solid angle) 



Radiated Power P rar j : 



Prad := 2 -ti- 



U(0)sin(0)d9 



(W) 



Prad = 1.256 



(W) 



Directivity D 0 : 

Umax •= 1 (W/solid angle) 



Do := 



4 n Umax 



Prad 



(dimensionless ) 



Do = 10.007 



(dimensionless ) 



Effective Isotropic Radiated Power EIRP: 



EIRP := Prad Do 



(W) 



EIRP = 12.566 



(W) 
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Radiation Resistance R 



r_i 



Rr : = 



2 Prad 

( |lo| ) 2 



( 0 ) 



Rr = 2.512 



(Q) 



Input Resistance R lrt : 

Rin := Rr (Q) 

Rin =2.512 (Q) 



Voltage Relection Coeffecient f: 

Rin - Zo 
Rin + Zo 



(dimensionless ) 



f =-0.935 



(dimensionless) 



Reflection Efficiency : 

et = ccd [ 1 - ( |f| ) 2 ] (dimensionless ) 



et =0 125 



(dimensionless ) 



Ga in G i 



G- = €t Do (dimensionless) GdB •= 101og(G) (dB) 



G = 1.255 (dimensionless) GdB =0.986 (dB) 



Maximum Effective Aperture : 



Aem := 



G-r 

4*71 



(m 2 ) 



Aem = 0. 1 



(m 2 ) 



Maximum Effective Height h ^ m : 



hem := 2- 



f Rr-Aem 

rio 



(m) 



hem = 0.052 



(m) 
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NORMALIZED ARRAY FACTOR RECTANGULAR PLOT 



Number of radiating elements: 



N = 10 



Direction of Main Lobe (radians): 
0o = 0 




0 0.628 1.257 1.885 2 513 3.142 3.77 4.398 5.027 5.655 6.283 

0 
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NORMALIZED ARRAY FACTOR POLAR PLOT 



Number of Elements: N - 10 

Direction of Main Beam (radians) : 0o = 0 



X(9):« | AFn(6)| • cos] 0 + “ 



Y(0) := |AFn(0)| •sin|0 + — 



0 = 0 ° 




0 = 180 ° 
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PLANAR ARRAY 

MATHCAD SOFTWARE ARRAY_PL . MCD 

This planar array application is based on M and N isotropic point 
sources uniformly spaced and aligned on the x-axis and y-axis 
respectively. Mutual interference between adjacent elements are 
neglected when calculating antenna parameters. Since isotropic point 
sources are used as the radiating elements, polarization of the antenna 
can not be determined. The gain can only be idealized when it is 
assumed that the antenna input resistance is equivalent to the radiation 
resistance and the total antenna efficiency is calculated for a ideal 
lossless antenna. 

(Note: Mathcad equations can not use symbolic subscripts. Therefore, 

parameters do not contain subscripts (i.e., AF n will be written AFn) . ) 



parameters (Items with * indicate parameters that can be 


computed but 


their 


results are idealized as the result of using isotropi 


c point 


sources as the radiating elements) : 








\ = 


Wavelength 








Lx = 


Antenna Dimension in x-direction 








L y = 


Antenna Dimension in y-direction 








r min 


= Minimum Distance to the Far-Field 








Px = 


Progressive Phase Shift Between Adjacent 


Elements 


in 


x-direction 


py = 


Progressive Phase Shift Between Adjacent 


Elements 


in 


y-direction 


Yx = 


Array Factor Phase Shift in x-direction 








Yy = 


Array Factor Phase Shift in x-direction 








AF n ( 0 , <j) ) = Normalized Array Factor 








U 


= Radiation Intensity 








p rad 


= Radiated Power 








a 

0 

ii 


Directivity 








EIRP 


= Effective Isotropic Radiated Power 








R r = 


Radiation Resistance 








R in 


= Antenna Input Resistance * 








r = 


Voltage Reflection Coefficient * 








e t = 


Antenna Efficiency * 








G = 


Gain * 








A em 


= Maximum Effective Aperture * 








h em 


= Maximum Effective Height * 
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The following known or estimated data must be entered: 



M = Number of Isotropic Radiating Elements in x-direction 
N = Number of Isotropic Radiating Elements in y-direction 
f = Frequency of Interest 

d x = Distance Between Adjacent Elements in x-direction 

dy = Distance Between Adjacent Elements in y-direction 

0 O = Direction of Main Beam at 0 = 0o 

<(> 0 = Direction of Main Beam at = tyo 

I Q = Antenna Feed Current 

Z Q = Characteristic Feed Impedence 

j = Number of Increments in Degrees for Far-Field Radiation Pattern 



M := 5 



Enter input data here 



( elements ) 



dx :=0.5 



(meters ) 



N =5 



(elements ) 



dy =0.5 (meters) 



f := 3 10 8 (Hz) 




(radians) 



Io =1.0 (A) 




( radians ) 



Zo:=75.0 (Q) 



j := 180 



(degrees ) 
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Define constants and calculate wavelength: 



c := 2.9979- 10 8 (meters /sec) 



c 

1 := - (meters/cycle) 

f 



1 =0.999 (meters/cycle) 



r|o := 1 20 n 




(Q) 



(m 1 ) 



ccd := 1 



(dimensionless ) 



Calculate linear array parameters: 
Define angular offset 0_ and jfe: 

e ; (o~ io' 4 ) io' 4 + — j . (2 n- 1CT 4 ) 

§ := (0- lOT^ ia 4 + — J .. (2 71 - icr 4 ) 



( radians ) 



( radians ) 



Minimum Distance to Far-Field r m i n _^ 
r Q := 1.61 (meters ) 

r, '--5 dx) 2 + (N dy) 2 (meters) 

2- 



u 



(M dx) 2 + (N dy) 2 



( m*=>t prcil 



nnin := max(r) 



(meters ) 



rmin = 25.018 



(meters ) 



Array Factor AF n : 

Px :=-k-dx-sin(6o)-cos(^o) ( radians) 



Py '=-k*dy sin(8o) sin((t>o) 



( radians ) 



yx(0,<t>) : = k-dx-sin(0)-cos(<|O-h Px 



( radians ) 



iyy(0,<|O ' = k-dy-sin(0)-sin(<(O + Py 



AFn(0,<(>) • = 



sinly VX(9,|) 

M-antep 




Radiation Intensity U: 



U(0,4O •= AFn(0»AFn(0 # 4) 



( radians ) 



(dimensionless ) 



(W/solid angle) 



Radiated Power P r<= 1 H : 
' 2-n 



Prad := 



0 



U(0,4>)sin(0) d0d<(> 



0 



(W) 



Prad =0.206 



(W) 
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(W/solid angle) 



Directivity D n : 

Umax := 1 



_ 4 n Umax 

Do .= 

Prad 



Do =61.09 



(dimensionless ) 



(dimensionless ) 



Effective Isotropic Radiated Power EIRP: 

EIRP = Prad Do (W) 

EIRP = 12.566 (W) 



Radiation Resistance R r : 

Rr: = i±^ 

( |lo| ) 2 



(Q> 



Rr =0.411 



(Q) 



Input Resistance R j n : 



Rin := Rr 



(Q> 



Rin =0.411 



(Q) 
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Voltage Relection Coeffecient Y: 

j, _ Rin - Zo 
Rin + Zo 

r =-0.989 

Reflection Efficiency £ t : 

et :=tcd [l - (|r|) 2 ] 

et =0.022 
Gain G: 

G: = ct Do (dimensionless) GdB := 101og(G) 

G = 1.326 (dimensionless) GdB = 1.225 

Maximum Effective Aperture 

r i 2 

Aem := (m 2 ) Aem =0.105 

4 71 

Maximum Effective Height h ^ m : 

hem := 2- (m) hem = 0,021 

4 i i° 



(dimensionless ) 



(dimensionless ) 



(dimensionless ) 



(dimensionless ) 



(dB) 



(dB) 



(m 2 ) 
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TWO-DIMENSIONAL ANTENNA PATTERN: (x-z) plane 



Number of 


Elements 


along x-axis: 


2 : 

H 




Number of 


Elements 


along y-axis: 


N =5 




Direction 


of Main 


Beam: 


00=0.524 


( radians ) 








= 0.785 


(radians) 



X1(0) •= |AFn(0,O)| cosjo- Y1(0) '= | AFn(0,O)| sin|o- ^ 

0 = 0 ° 
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TWO-DIMENSIONAL ANTENNA PATTERN: (y-z) plane 



Number of 
Number of 
Direction 



Elements along 
Elements along 
of Main Beam: 



x-axis : 
y-axis : 



X2(0) 



AFn 




•cos 0- 



M =5 
N =5 
Oo =0.524 
<t>o =0.785 

Y2(0) = 



( radians ) 
( radians ) 




e = o° 
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TWO-DIMENSIONAL ANTENNA PATTERN: $ = fl/4 



Number of 
Number of 
Direction 



Elements along 
Elements along 
of Main Beam: 



x-axis : 
y-axis : 



X3(0) 



AFn 




M =5 
N =5 
Go = 0.524 
(K> = 0.785 

Y3(0) : = 



( radians ) 
( radians ) 




0 = 0 ° 
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TWO-DIMENSIONAL ANTENNA PATTERN: <fe = $ Q 



Number 


of 


elements 


along x-axis : 


M =5 


Number 


of 


Elements 


along y-axis: 


N =5 



Direction of Main Beam: 00 = 0.524 (radians) 



4>o = 0.785 (radians) 

X4(0) := |AFn(0,$o)|-cos|e--iJ Y4(0) := | AFn(0,*o)| sin^0- 

0 = 0 ° 




xref ,X4( 6) 
m 
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CIRCULAR ARRAY 

MATH CAD SOFTWARE ARRAY CI.MCD 



This circular array application is based on N isotropic point sources 
uniformly spaced on the x-y plane along a circular ring of constant 
radius. Mutual interference between adjacent elements are neglected 
when calculating antenna parameters. Since isotropic point sources are 
used as the radiating elements, polarization of the antenna can not be 
determined. The gain can only be idealized when it is assumed that the 
antenna input resistance is equivalent to the radiation resistance and 
the total antenna efficiency is calculated for a ideal lossless antenna. 



(Note: Mathcad equations can not use symbolic subscripts. Therefore, 

parameters do not contain subscripts (i.e., AF n will be written AFn) . ) 



The circular array antenna Mathcad application will compute thefollowing 
parameters (Items with * indicate parameters that can becomputed but 
their results are idealized as the result of usingisotropic point 
sources as the radiating elements) : 

X = Wavelength 

r min = Minimum Distance to the Far-Field 

a n = Phase Excitation of the nth Element 

AF (0, <(>) = Normalized Array Factor 
U (0,<t>) = Radiation Intensity 

p rad = Radiated Power 

D q = Directivity 

EIRP = Effective Isotropic Radiated Power 
R r = Radiation Resistance 

R in = Antenna Input Resistance * 

T = Voltage Reflection Coefficient * 

e t = Antenna Efficiency * 

G = Gain * 

A em = Maximum Effective Aperture * 
h em = Maximum Effective Height * 
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The following known or estimated data must be entered: 



N = Number of Isotropic Radiating Elements 
f = Frequency of Interest 
a = Radius of circle 

0 O = Direction of Main Beam at 0 = 0o 

<)> 0 = Direction of Main Beam at $ = 

I Q = Antenna Feed Current 

Z Q - Characteristic Feed Impedence 

j = Number of Increments in Degrees for Far-Field Radiation Pattern 



Enter input data here 



N := 10 



(elements ) 



_10 
2 71 



(meters ) 



f:=3 10 8 



(Hz) 



Go := 0 



( radians ) 



Io := 1.0 



(A) 



<(>o := 0 



( radians ) 



Zo :=75.0 



(Q) 



j •= 360 



(degrees ) 



Define constants and calculate wavelength: 



c := 2.9979- 10 8 



(meters/sec) 



t|o := 1 20 n 



(Q) 




f 



(meters/ cycle) 
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Define constants and calculate wavelength: 



31 =0.999 (meters/cycle) ecd := 1 

Calculate linear array parameters: 

Define angular offset 0_ and _$: 

9 := (o- 10r 4 ) , 10 ‘ 4 + — J„ (2 7 i- icr 4 ) 
t = (0- icr 4 ),|- io‘ 4 + (2 7 t- icr 4 ) 



Minimum Distance to Far-Field r m i n : 



r Q = 1.6 X (meters) = 5(2a) 



2 (2 a) 2 
2 " X 



(meters ) 



rmin := max(r) (meters) 



rmin = 20.278 



Array Factor AF: 
n = 1 .. N (elements) 



(I) = 2 71 — 

n N 



: -k a sin(0o) cos^<j>o- d)^ 



(dimensionless ) 



( radians ) 



( radians ) 



(meters ) 



(meters ) 



( radians ) 



( radians ) 
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Array Factor AF : 



AFce.f) 




■ a s»n( 0)-cos(<t>- ^+0^ 



n 



(dimensionless ) 



Radiation Intensity U: 



U( 0, <)>) := AF(0,*)AF(0,*) 



(W/solid angle) 



Radiated Power P rar ) : 



2*71 



71 

2 



Prad := 



♦0 



U(G,4>) sin(0) d0 d<)> 

• 0 



(W) 



Prad =0.536 



(W) 



Directivity D n : 



Umax : = 1 



(W/solid angle) 



Do •= 



4- 7i- Umax 
Prad 



( dimensionless ) 



Do =23.427 



(dimensionless ) 



Effective Isotropic Radiated Power EIRP: 



EIRP := Prad Do 



(W) 



EIRP = 12.566 



(W) 
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Radiation Resistance R r : 

„ 2-Prad 
Rr := 

(|lo| ) 2 



<fl> 



Rr= 1.073 



(Q) 



Input Resistance R j n : 

Rin := Rr (Q) 

Rin = 1.073 (Q) 



Voltage Relection Coeffecient f: 

P Rin - Zo 
Rin + Zo 



(dimensionless ) 



r =-0.972 



(dimensionless) 



Reflection Efficiency : 

et = ecd [ 1 - ( lrl) ; ] (dimensionless) 



et =0.056 



(dimensionless) 



Gain G : 



G=et-Do (dimensionless) GdB — 101og(G) (dB) 



G= 1.303 (dimensionless) GdB = 1.149 (dB) 



Maximum Effective Aperture 

G . i 2 

Aem := (m 2 ) Aem =0.104 (m 2 ) 

4*71 



Maximum Effective Height h ^ m : 

hem =2- j fi 1 ^ cm ( m ) hem =0.034 (m) 

N VP 
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TWO-DIMENSIONAL ANTENNA PATTERN: (x-z) plane 



Number of Elements in x-y plane: 
Direction of Main Beam: 



X 1 ( 6) := |AF(e,0)| •cos|g- 



N = 10 

0o=O (radians) 

<(>o=0 (radians) 

Y 1(0) := |AF(0,O)| • sin^0 — — J 

e = o° 






Y1(0) 



xref ,X1(9) 
m 




126 



TWO-DIMENSIONAL ANTENNA PATTERN: (y-z) plane 



Number of Elements in x-y plane: 
Direction of Main Beam: 



X2(0) 




N = 10 
0o=O 
4>o =o 

Y2(0) 



( radians ) 
( radians ) 




0 = 0 ° 




xref m , X2{ 0) 



127 



TWO-DIMENSIONAL ANTENNA PATTERN: & = n / 4 



Number of Elements in x-y plane: 
Direction of Main Beam: 



X3(0) 




N = 10 

0o=O (radians) 

4>o = 0 (radians) 



Y3(0) 




0 = 0 ° 
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TWO-DIMENSIONAL ANTENNA PATTERN: & = 



Number of Elements in x-y plane: 
Direction of Main Beam: 



X4(0) := |AF(0,*o)| cosfe- i 



N = 10 




0o=O 


( radians ) 


4>o = o 


( radians ) 



Y4(0) := | AFC0,*o)| - sinf 0 — - 



0 = 0 ° 
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FOLDED DIPOLE 

MATHCAD SOFTWARE FOLDED. MOD 



To provide impedance matching to coaxial input cables and maintain good 
directional pattern characteristics, a folded dipole is used instead of 
a single dipole element. The geometry of a folded dipole is such that 
the spacing between dipole elements (d) is much less than the wavelength 
(1) and the length of an element (L) . 

This Mathcad application calculates the parameters of a folded dipole 
aligned vertically on the z-axis and a folded dipole positioned 
horizontally above the earth. For the calculations , the antenna is 
assumed to behave as a cylindrical dipole with equivalent radius (a e ) . 

(Note: Mathcad equations can not use symbolic subscripts. Therefore, 

parameters do not contain subscripts (i.e., Rj^ n will be written Rin) . ) 

The folded dipole antenna Mathcad application will compute the following 
parameters (* indicates antenna parameters that are calculated for a 
folded dipole in free space and positioned horizontally over the earth) : 



X = Wavelength 

a e = Equivalent Radius 

R r = Radiation Resistance 

R^ = Input Resistance of Linear Dipole 

X r = Reactance of Linear Dipole 

= Input Reactance of Linear Dipole 
Z^ = Input Impedance of Linear Dipole 

Z Q = Characteristic Impedance of Two-wire Transmission Line 
Z t = Characteristic Impedance Transform of Folded Dipole 
Z^ n = Characteristic Impedance of Folded Dipole 
Z Q = Maximum Current of Folded Dipole 

r min = Minimum Distance to the Far-Field 
E { 0 ) - Electric Field * 

E n {0) = Normalized Electric Field * 

U(0) = Radiation Intensity * 

Umax = Maximum Radiation Intensity * 

P rac j = Radiated Power * 

D 0 = Directivity * 

EIRP = Effective Isotropic Radiated Power * 

Rin = Antenna Input Resistance * 

f = Voltage Reflection Coefficient * 
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8 t = Antenna Efficiency * 

o a = Antenna Polarization Vector * 

PLF = Polarization Loss Factor * 

G = Gain * 

A em = Maximum Effective Aperture * 
h em = Maximum Effective Height * 

s r p = Relative Complex Effective Dielectric Constant of 
T v = Vertical Reflection Coefficient of Ground Plane 

= Horizontal Reflection Coefficient of Ground Plane 



Ground 



The following known or estimated data must be entered: 

f = Frequency of Interest 

a = Radius of Conductor 

b = Radius of Conductor 

d = Distance Between Conductors 

L = Length of Antenna 

h = Height of Antenna Above Earth 

e cc j = Conduction/Dielectric Efficiency 

e r = Relative Dielectic Constant of Ground 

0 = Conductivity of Ground 

r ff = Far Field Distance 

1 = Number of Increments in Degrees for Far Field Radiation Pattern 
o w = Incoming Wave Electric Field Unit Vector For Free Space Antenna 



°wh 

0 ph 



INcoming Wave Electric Field Unit Vector For Horizontal Antenna 
Antenna Polarization Direction For Horizontal Antenna 



Enter input data here 





i 






'0' 


aw := 


0 


(dimensionless ) 


owh •= 


0 




loj 






ll. 



(dimensionless ) 



f := 3 10* (Hz) 



0ph :=0 



( radians ) 



a := 0.0005 (meters ) 



rff := 1 ■ 10 5 



(meters ) 



b ’= a 



(meters ) 



i := 90 



( radians ) 



d := 12. 5 a (meters ) 



ecd := 1 



(dimensionless ) 



L := 0 5 



(meters ) 



o:=25 



( siemens /meter ) 



h := .5 



(meters ) 



a* := 15 



(dimensionless) 
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Define constants: 



a a 



V 

o 

w 



(dimensionless ) 



c := 2.9979- 10 8 (meters/sec) 



TjO := 1 20 71 



X 



c 

f 



(meters /cycle) 




X = 0.999 



(meters /cycle) 



7 := 0.57721 



£0 • 



— iff 9 
36 n 



(F/ra) 



z- = 0..2-i 



Define angular offset 8: 



i:=0.— ..2-Jt 

i 



Define cosine [Ci(x)] and sine [Si(x)] integrals: 



n := 1 .. 50 



Ci(x) =0.57721 +• 



ln(x)~y^ 



(- 1)"" 1 x 2n 
(2 n) (2 n)! 



n 



Si(x) -=2 



(-D 



(2 n- 



n 



(fl) 



(m !) 



(dimensionless ) 



( radians ) 



(radians ) 



n- 1 x < 2- n — 1) 

1 ) (2 n- 1)! 
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FOLDED DIPOLE ALIGNED ON THE Z-AXIS IN FREE SPACE: 



Equivalent radius a ^: 



ae - exp 



a 2 ln(a) + b 2 ln(b) + 2 a b ln(d) 



(a+ b) 



ae = 0.002 



Electric field E ( 9 ) : 



E( 0 ) =j T)o * 



-j k rff 
e J 

2 Ti rff 



^cos(-^ cos(0)| - cost^ 

l 2 } \2l 



sin( 0 ) 



En( 0 ) := 1 ^- 



: (l 



Minimum Distance to Far-Field 



mm_L 



r 0 :=l 6 X 



(meters ) 



r, :=5 L 






(meters ) 



rmin:=max(r) (meters) 



nmm = 2.5 



(meters ) 



(meters ) 



(meters ) 



(V/m) 



(V/m) 



(meters ) 
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Radiation Intensity U ( 0 ) : 



U(0) 



■ ( |E(6)| ) 2 

2 TIO 



(W/solid angle) 



Radiated Power P r ^ H : 

"71 

Prad := 2 -ti- U(0) sin(0) d0 

• 0 



(W) 



Prad = 36.64 



(W) 



Directivity D n : 



Umax :=U 




(W/solid angle) 



„ 4-7i Umax 

Do := 

Prad 



(dimensionless ) 



Do = 1.641 



(dimensionless ) 



Effective Isotropic Radiated Power EIRP: 



EIRP := Prad- Do 



(W) 



EIRP =60.132 



(W) 
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Radiation Resistance R r : 



Rr : = — 
2n 



y-f ln(k L) - Ci(k L) + — sin(k L) (Si(2 k L) - 2 Si(k L)) 
2 



+ — cos(k-L) 
2 



HtI 



4- Ci( 2 k L) - 2 Ci(kL) 



(Q) 



Rr = 73.281 



(Q) 



Input Radiation Resistance of linear dipole : 




Rd =73.281 (Q) 



Reactance of linear dipole X r : 



Xr: = ^ 
4 n 



2 Si(k L)+ cosfk L) (2 Si(k L)- Si(2 k L)) 

1 sin(kL) 2 Ci(kL)- Ci(2 kL)- Ci p k ae 



1 ) 



(Q) 



Xr = 43.142 



(Q) 



Input reactance of linear dipole X^-i : 




Xd =43.142 (Q) 
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Input impedance of linear dipole Z^: 



Zd = Rd+j Xd 



(Q) 



Zd =73.281 +43.142i 



(£ 1 ) 



Characteristic impedance of two-wire transmission line Z n : 
Zo: = — acosh] — - — j (Q) 

* \2 JTb/ 



Zo =302.312 



(Q) 



Characteristic impedance transform for folded dipole Z t : 



Zt - j 




(D) 



Zt =-2.747 10 5 i 



(Q) 



Characteristic impedance of folded dipole Z -j n : 

. 4 Zt Zd 

Zin . = 

2Zd+ Zt 



Zin =293.307 +172.465i 



(Q) 



Input Resistance R j n : 



(Q) 



Rin =Re(Zin) 



( 0 ) 



Rin =293.307 



Voltage Relection Coeffecient £ : 



_ Rin - Zo 
Rin + Zo 


(dimensionless) 


r = -o.oi5 


( dimensionless ) 


Reflection Efficiency : 






«t = ecd [ 1 - (|r|) 2 ] 


(dimensionless) 


et = 1 


( dimensionless) 


Gain G i 








G: = et Do 


( dimensionless ) 


GdB := 101og(G) 


(dB) 


G = 1.641 


(dimensionless ) 


GdB = 2.15 


(dB) 


Polarization loss factor PLF : 






PLF := flow oa|) 


(dimensionless ) 


PLF = 1 


(dimensionless ) 



Maximum Effective Aperture 



A • T> 

Aem •= PLF 


- -erru. 

(m 2 ) 


Aem = 0. 1 3 


(m 2 


4*71 








Maximum Effective Height h^ m : 






hem -2* K Aem 


(m) 


hem =0.318 


(m) 


4 T]0 
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Far Field radiation pattern of vertical folded dipole in free space. 



X(0) := |En(0)|-cos|O+-2 



Y(0):= |En(0)| sinjo+^j 



0 = 0 ° 




0 = 180° 
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HORIZONTALLY POSITIONED FOLDED DIPOLE OVER THE EARTH: 



Define angular offsets 0 1 and 02: 



01 + JL..JL (radians) 

2 2 2i 2 



02 :=-i + JL, 

z 2 2-i 



Relative complex effective dielectric constant g r p j_ 



crp := er- j 



2-n-feo 



^ component of E-field E j>: 



E(K01 ) : = j rp- 



-1 k iff 

e 3 

rvA. 


cos|^-^ sin(Gl )sin(0)J - cos|— J 


rp - 

2-n-rff 


r ' I 2 




-^1 - sin(01 ) sm(0) 2 



Eijil :=j ■ t]o~ 



e-j krfr 


cos|-^ sin^G2^ sin(0)j - cos|^-ij 


2-jirff 


- sin^G2^ 2 sin(O) 2 



6 component of E-field EQ: 




( radians ) 



(dimensionless ) 



(V/m) 



(V/m) 



(V/m) 



(V/m) 
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Vertical plane wave reflection coefficient f w : 



rv(01 ) .. ap cos(ei)- a/oT- sin(Ql ) 2 

erpcos(01 )4-^Jerp- sin(01) 2 



fop- sin ^02^ 2 

opcos(02^ 4 J&p- sin ^02^ 2 



erp* cos| 



N-J 



(dimensionless ) 



(dimensionless ) 



Horizontal plane wave reflection coefficient £ h » 



j cos(91 )- -Jap - sin(6I) 2 

( 2 

cos(01 ) 4 ^erp- sin( 0 1 ) 



(dimensionless ) 



cos (02 \ - Jerp - sin (02 ) 2 

rhl * = - — 2i 1 — _ (dimensionless) 

cos ^02^ 4 ^fap - sin ^02^ 2 



E-field ( ^ component) with ground plane reflection E j> gp J_ 
E4>gp(01 ) •= E4>( 91 )• (e* k h cos(ei) + rh(01) e* j k h cos( 61 >) (V/m) 



E«(>gp 1 2 : = 



[e 4> 1 z (e 



k-h-cos 



H 



rhl e 

z 



-j kh-cos 



H) : 



(V/m) 



max(E<(igpl) 



(V/m) 
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E-freld ( 6 component) with ground plane reflection EQgp 



E0gp(01 ) :=E9(01 )• (e J k h cos(6l) - fv(01 ) e‘ j k h ' cos(6l) ) 



(V/m) 



E0gpl z : = 



E01 



(e* k h ' 



■(») 



- fvl e 



khcos 



to) 



(V/m) 



E0n(01 ) := ^ R9gp( ^ 1 
max(E6gpl ) 



(V/m) 



Radiation Intensity Uh( 6 l): 



o 

Uh( 0 1 ) =— [( |E<frgp(0l )| ) 2 + ( I E0gp( 0 1 )| ) 2 ] 

2rp> 



(W/solid angle) 



Radiated Power Ph rar j: 



Phrad 2-7i 



Uh(01 )sin(01 ) d01 



(W) 



Phrad = 70.415 



(W) 



Directivity Dh n : 



Ugp_max z 



iff 2 

2-rp 





(W/solid angle) 



Uhmax = max( U gp _max ) 



(W/solid angle) 
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Directivity Dh 0 : 



4-7i Uhmax . , . , 

Dho = (dimensionless) 

Phrad 



Dho = 3.938 (dimensionless ) 



Effective Isotropic Radiated Power EIRPh: 

EIRPh := Phrad Dho ( W ) 



EIRPh =277.3 



(W) 



Input 
Rhin := 



Resistance 



2 Phrad 




Rh 



m _ l 



Rhin = 140.83 



(Q) 



(Q) 



Voltage Rejection Coeffecient T H : 

fH := ^ in ~ (dimensionless ) 

Rhin Zo 



fH =-0.364 (dimensionless) 



Reflection Efficiency eh h : 

eht : = ecd [l - ( |fH| ) 2 ] (dimensionless) eht =0.867 (dimensionless) 
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Gain Gh : 



Gh -eht-Dho (dimensionless) GhdB - 10 log(Gh) (dB) 



Gh =3.415 (dimensionless) GhdB =5.334 (dB) 



Polarization loss factor PLFv: 
Ehmax’ = | E0n( Gph) | (V/m) 



Lobes ‘= floor 




aah : = 



0 

0 

Ehmax 

Lobes 



(dimensionless ) 



PLFh : 



.(i 



owh oah|) 



(dimensionless ) 



PLFh =7.736 10 



(dimensionless ) 



Maximum Effective Aperture Av ^ m : 
r'h i ^ 

Ahem — PLFh (m^) Ahem = 2.099" 10 4 (m^ ) 

4-71 



Maximum Effective Height hv ^: 

,, . . [Rhin Ahem 

hhem • = 2 / — (m) 



tio 



hhem =0.018 (m) 
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H-Plane Radiation Pattern (y-z plane) of Horizontal Folded Dipole Over 
The Earth: 



Y 1 ( 01 ) = | E4>n( 0 1 )| 



sin 01 




Xl(01 ) = | E«)>n( 0 1 )| 



cos 1 01 + - 




xref m ,Xl( 01 ) 
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E-Plane Radiation Pattern (x-z plane) of Horizontal Folded Dipole Over 
The Earth: 



Y2(01 ) - | E0n( 0 1 )| •sin 




X2(01 ) '= |E0n(01)| 





xref , X2( 01 ) 
m 
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CAGED DIPOLE ANTENNA 
MATHCAD SOFTWARE CAGED_DI.MCD 

This caged dipole application is modeled after a circular array with the 
radiating elements as dipoles. The dipoles are modeled after conductors 
with infinitely small radii and antenna length greater than 1/4. Mutual 
interference between adjacent elements are neglected when calculating 
antenna parameters. 

(Note: Mathcad equations can not use symbolic subscripts. Therefore, 

parameters do not contain subscripts (i.e., AF n will be written AFn) . ) 



The caged dipole antenna Mathcad application will compute the following 
parameters : 

VERTICAL CAGED DIPOLE IN FREE SPACE: 

1 = Wavelength 

r min ” Minimum Distance to the Far-Field 

AF (G^) = Normalized Array Factor 
E (0,<t>) = E-Field of Finite Length Dipole 

E t (0,<|)) = Total E- Field 

E tn (0i) = Total Normalized E-Field 

U (0,<|>) = Radiation Intensity 

P ra d = Radiated Power 

D q = Directivity 

EIRP = Effective Isotropic Radiated Power 

R r = Radiation Resistance 

R in = Input Resistance 
F = Voltage Reflection Coefficient 
e t = Antenna Effeciency 
G — Ga i n 

PLF = Polarization Loss Factor 
A em = Maximum Effective Aperture 

hem = Maximum Effective Height 
flow = Low Frequency 

fhigh = High Frequency 
BW = Bandwidth 

VERTICAL CAGED DIPOLE OVER EARTH: 

8 rp = Relative Complex Effective Dielectric Constant 
AFv(01,<|O = Normalized Array Factor 
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Ev (01 , <j>) = E-Field of Finite Length Dipole 

T v (81 ) = Vertical Plane Wave Reflection Coefficient 

Ev t (01 , $) = Total E-Field 

Ev tn (01,<(>) = Total Normalized E-Field 

Uv(01,<t») = Radiation Intensity 

Pv rad = Radiated Power 
Dv 0 = Directivity 

EIRPv = Effective Isotropic Radiated Power 
Rv r = Radiation Resistance 

R v in = Input Resistance 
rv = Voltage Reflection Coefficient 
ev t = Antenna Efficiency 
Gv = Gain 

PLFv = Polarization Loss Factor 
Av eiT1 = Maximum Effective Aperture 

hv em = Maximum Effective Height 

HORIZONTAL CAGED DIPOLE OVER EARTH: 

AFh (01 , <(il ) = Normalized Array Factor 

1^(01) = Horizontal Plane Wave Reflection Coefficient 

E0h(01) = 0 Component of E-Field 

E0h t (01 , <j»l ) = 0 Component of Total E-Field 

E0h tn (01, ()il ) = 0 Component of Total Normalized E-Field 

E<|>h(01) = <J» Component of E-Field 

E<t»h t (01, $1 ) = $ Component of Total E-Field 

E^hfn (8 ^ t H ) = Component of Total Normalized E-Field 

Uh (01 , <J>1 ) = Radiation Intensity 

ph rad = Radiated Power 
Dh 0 = Directivity 

EIRPh = Effective Isotropic Radiated Power 
Rh r = Radiation Resistance 

Rh in = Input Resistance 

fH = Voltage Reflection Coefficient 
eh t = Antenna Efficiency 
Gh = Gain 

PLFh = Polarization Loss Factor 
Ah e m = Maximum Effective Aperture 

hh eltl = Maximum Effective Height 
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The following known or estimated data must be entered: 



N = Number of Dipoles 

a = Radius of Antenna 

L = Length of Antenna 

h = Antenna Height Above Ground 
f = Frequency of Interest 
I Q = Antenna Feed Current 

Z Q = Input Impedance 

e r = Relative Dielectric Constant of Ground 
g = Conductivity of Ground 

a w = Incoming Wave E-Field Unit Vector for Antenna (Free Space) 

o a = Antenna Unit Polarization Vector (Free Space) 

ov w = Incoming Wave E-Field Unit Vector for Vertical Antenna 

0Vp = Direction of Incoming Unit Vector (Vertical Antenna) 

gh w = Incoming Wave E-Field Unit Vector for Horizontal Antenna 

9hp = Direction of Incoming Unit Vector (Horizontal Antenna) 

j = Number of Increments in Degrees for Far Field Radiation Pattern 



Enter input data here: 





IA 






M 


gw := 


0 


(dimensionless ) 


ga := 


0 




lol 






,0, 



(dimensionless ) 



ovw - 



(dimensionless ) 



_ 80n 

0vp = 

180 



(dimensionless ) 



ahw := 



(dimensionless) 0hp:=O (dimensionless) 



N =8 



(elements ) 



Io = 1.0 



(A) 



a := - 

8 



Zo : = 75.0 



(meters ) 
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(Q) 



L := 0.5 



(meters ) 



j := 90 



(increments ) 



h •= 0.25 



(meters ) 



c •= 10 2 



f := 3 10 8 (Hz) 



a- := 15 



Define constants and calculate wavelength: 
c := 2.9979- 10 8 (meters/sec) r |0 := 1 20* ti 



k 



c 

f 




to := — — 10' 9 
36 n 



(meters/cycle ) 



(A) 



( farads/m) 




rff ■- 10 5 



ccd := 1 



k = 0.999 



(meters/cycle ) 



z •= 0 .. 2 j 



Calculate caged dipole parameters parameters: 
Define angular offset 6_: 



e =0,— .2 n 
j 



( radians ) 



4’rO.— ,.2-JI 

j 



( siemens/m) 



(dimensionless ) 



(fi) 



(m 1 ) 



(meters ) 



(dimensionless ) 



(increments ) 



( radians ) 
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Minimum Distance to Far-Field r m in : 



r Q := 1.6 X (meters) 



r,:=5^ 



L 2 4 (2 a) 2 



(meters ) 



2-[l 2 4 (2 a) 2 ] 

l 



(meters ) 



rmin := max(r) (meters ) 



rmin = 2.795 



VERTICAL CAGED DIPOLE IN FREE SPACE: 



Array Factor AF ( 6 ) : 

n •= 1 .. N (elements) :=2 n — 

n N 



AF(O^) :=y^In e 



kasirH 0) cos 



E-Field of finite dipole E ( 6 ) : 



E(6) :=J 



Tp 



T -i k rflf 
In-e J 


cos| 


^. CO S(0)' 


| - cos| 


m 


2 Tt rflf 


i sin(0) 



(meters ) 



( radians ) 



(dimensionless ) 



(V/m) 
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Total E-Field E j. ( 6) & Normalized E-Field E m (6) : 



Et(0,(t>) := E(0) AF( 9, 4>) 



(V/m) 



Em(9, <f>) 



Eue,^) 




Radiation Intensity U( 9 ): 

U(0,<t») •=——•( | Et( 0, (t>)| ) 2 
2t]0 



(V/m) 



(W/solid angle) 



Radiated Power P 



' 2 • n 71 



Prad := 



•0 



U(0,<j>) sin(0) d0d<j> 

•0 



(W) 



Prad =0.441 



(W) 



Directivity : 



Umax : = U 




(W/solid angle) 



4 7i- Umax 
Do := 

Prad 



(dimensionless ) 



Do = 1.543 



(dimensionless ) 



Effective Isotropic Radiated Power EIRP: 



EIRP := Prad Do (W) 



EIRP =0.681 (W) 
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Radiation Resistance R r : 

0 • PrflH 

Rr :z—LL— (Q) Rr = 56.487 (Q) 

(|ln |) 2 



Input Resistance R in : 



Rin := (Q) Rin = 56.487 (Q) 




Voltage Relection Coeffecient f: 

F:_Rin — Zo (dimensionless) T =-0 141 (dimensionless) 

Rin + Zo 

Reflection Efficiency g t : 

Et '= Ecd [ 1 - (|r|) 2 ] (dimensionless ) et=0.98 (dimensionless ) 

Ga in G ! 

G- = et-Do (dimensionless) GdB •= lO-log(G) (dB) 

G = 1.513 (dimensionless) GdB = 1.798 (dB) 



Polarization Loss Factor PLF : 






PLF 



(dimensionless) 



PLF = 1 



(dimensionless ) 
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Maximum Effective Aperture 
Cl 2 

Aem := PLF (m 2 ) Aem=0.12 (m 2 ) 

4-71 



Maximum Effective Height : 

hem := 2 j ^ ^ em ( m ) hem =0.268 (m) 

4 T]0 



Define Frequency Relationship F(x): 

F(x) = - 5.696 10’ 5 x + 0.3 1481 (dimensionless) 



Frequency Range f low and f high j_ 



flow := f- 




•f 



(Hz) 



flow = 2.056 10 8 ( Hz ) 



fhigh :=f+ F 




f 



(Hz) 



thigh =3 944 1 0 8 (Hz) 



Bandwidth BW: 

BW = thigh - flow (Hz) BW = 1 .888* 10* (Hz) 
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E -FIELD PATTERN 



: This far-field 

caged dipole in free space: 

Number of Dipole Conductors: 

Radius of Antenna: 

X(0).= |Etn(e,*)| •cos|0- —I 



radiation pattern is for the vertical 
N =8 

a =0.125 (meters) 

Y(9) •= |Etn(0,*)| - sin^G — — | 

0 = 0 ° 




0 = 180 ° 
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VERTICAL CAGED DIPOLE OVER EARTH: 



Define angular offset 9 1 and _0 2 : 

( radians ) 



01 =-1,-1+ 71 71 



2 2 2 j 2 



4>2 - — • z 

z j 



62 =-i4JL. z 
2 2 2 j 



( radians ) 



Relative complex effective dielectric constant c r p ^ 
o 

crp = er- j - 



2-n-f-eo 

Array Factor AFv ( 6 ) : 

S i ka-sin(0I ) cos(ij>_ Q.) 

Ine V * 

n 

Afvgp i 

n 

E-Field of finite dipole Ev( 0 1): 



Ev(01) =j tp 



In e 



j k-rff 



2-ti rff 



. k L ... A /k-L\\ 

cos| cos(01)l - cosl — 



sin( 01 ) 



Evgp, : =j no 



, -j k-rff 

In e J 
2 7i rff 



cos 



fr co i >2 )) - 



sin^02^ 



Vertical plane wave reflection coefficient f v (0 1 ) : 



rv (0i ) - ap cos(ei ) - Jcrp- sin( 61 )‘ 



erp cos(01 )+■ -Jap- sin(01 )“ 



( radians ) 



(dimensionless ) 



(dimensionless ) 



(dimensionless ) 



(V/m) 



(V/m) 



( dimensionless ) 
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rvi 

z 



erp cos| 


M - J 


[op- sin | 


M 2 


erp cos) 


( e2 ) * J 


fop- sin | 


K ) 2 



Total E-Field Ev t (Q l) & Normalized E-Field Ev t- n (Q l): 
Evt(01 ,$) = Ev(01)AFv(01,i|>)(e i k h cos(6l> 4 Tv(01 ) e' j khcos(0l) ) 



Evtgp z 



Evgp z AFvgp z 




+• fvl 



-j k h cos 

•e 



H 



) 



E«ln(ei ,♦) = |Evl ! ol -*> | 

max(Evtgp) 



Radiation Intensity Uv( 6l) 

Uv( 0 U); s — ( |Evt(01 ,<t>)| ) 2 
2 TjO 



Uvl ( (Evlgpl) 

z 2 no 



Radiated Power Pv r ^d : 
’2-71 



Pvrad := 



0 



Uv(01 , (j))* sin( 01 ) dOl (fy 



0 



Pvrad = 0.483 



Directivity Dv 0 : 
Uvmax := Uv(0vp,O) 



(dimensionless ) 



(V/m) 



(V/m) 



(V/m) 



(W/solid angle) 



(W/solid angle) 



(W) 



(W) 
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(W/solid angle) 



Directivity Dv n : 
4 7i Uvmax 



Dvo := 



Pvrad 



Dvo = 4.4 15 



Effective Isotropic Radiated Power EIRPv: 

EIRPv := Pvrad Dvo (W) EIRPv =2.133 



Radiation Resistance Rv r : 

Rvt := 2 Pvr ‘ td (fi) Rvr =61.857 

(I Ini ) 2 



Input Resistance Rv .j n : 

Rvin '= — — — (fi) Rvin =61.857 




Voltage Rejection Coeffecient F V : 

rV := - ^ -- m — — (dimensionless) rV =-0.096 

Rvin 4- Zo 

Reflection Efficiency tv r : 

evt : = ccd [l - ( |fv| ) 2 ] ( dimensionless ) cvl =0.991 



(dimensionless ) 



( dimensionless ) 



(W) 



( 0 ) 



(Q) 



(dimensionless ) 



(dimensionless ) 
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Gain Gv : 



Gv =evtDvo 



(dimensionless ) 



GvdB = 10 log(Gv) (dB) 



Gv =4.374 



(dimensionless) GvdB =6.409 (dB) 



Polarization Loss Factor PLFv: 
Evmax := Evtn(0vp,0) ( V/m) 



ava ■ 



Evmax 

lobesv 

0 

0 



(dimensionless ) 



lobesv := floor 





PLFv 




(dimensionless ) 



PLFv =0.986 (dimensionless) 



Maximum Effective Aperture 

Gv • X. o 9 

Avem := PLFv (m 2 ) Avem= 0.343 (m 2 ) 

4 n 



Maximum Effective Height hv pm : 

hvem-= 2 - fc vr ^ve m ( m ) hvem = 0.474 (m) 

1 T]0 
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E-FIELD PATTERN: This far-field radiation pattern is for the vertical 

caged dipole over earth: 



Number 


of 


Dipole Conductors: 


N =8 


( elements ) 


Radius 


of 


Antenna : 


a =0.125 


(meters ) 


Height 


of 


Antenna Above Earth: 


h = 0.25 


(meters ) 



Xl(01) := |Evtn(01,<ti)|- cosOI + i 



Yl(01 ) := |Evtn(01,()o| •sin 




G = 0° 




xref , X 1 ( 01 ) 
m 
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HORIZONTAL CAGED DIPOLE OVER EARTH: 



Define angular offset $ 1 and J> 2 : 



41 *=0,-..2 n (radians) 

j 



Array Factor AFh( 9 l,^l): 
n *= 1 .. N ( elements ) 



0 '=2-71- 



N 



V 1 j k a fsin( <{>1 )-sin( 01 J-sinf©^ 4 -cos( 01 J-cosf©^) 

AFh(ei,*l)-=> lne V v * 

n 



AFhl z :=^In e 
n 



j k- 3 ' (sin ^$ 2^1 sin ^ 02^1 sin(©^ -+cos( 02 ^ cos(©^ 



9 component of E-Field of finite dipole EQh( Q l): 

fkl\ 



E0h(01 ) '= j i]0 



T -j krff 
lne 3 

2 7i rff 



/kL . 

)s sir 

\ 2 



cos| — - — sin( 0 1 ) | - cos| — 






sin(01 ) 



E0hl *=j • tjo 



In e’ 



j k rff 



2-ti rff 



cos| 


(T sin W 


J - cos| 


(t) 


J 


r — — 


Vn \ ^ 



Total E-Field E 9h f (01,$l): 

E0ht(ei,(|>l) :=E0h(01) AFKOl^l) (e 1 khcos < ei >_ rv(01) e' j khcos(01) ) 



EOhtl := 



E0hl AFhl 

Z Z 



| e J k h cos (®^ _ Fv ^ e ‘j ■ k 'h-cos(02^j | 



( radians ) 



(dimensionless ) 



(dimensionless ) 



(V/m) 



(V/m) 



(V/m) 



(V/m) 
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Normalized E-Field E 0h t n ( Ql , j) l ) : 

Eeiimiei ,n '■ |Eeh ‘ (el<t ' l > 1 

max(E0htl ) 



Horizontal plane wave reflection coefficient £^(01) 



cos(01 )- 


tap - sin( 01 ) 2 


cos(01 )+ aJ 


tap- sin(01 ) 2 



Thl : = 

z 



cos| 


W-J 


tap - sin| 




cos| 


KH 


ferp- sin| 


M 2 



$ component of E-Field of finite dipole E ^ h(8l) : 



In e J krfr 

E<t»h( 0 1 ) =J rp — | 1-cos 

2 Ji rff 



1 



In e"-* krff 

E4>h 1 '■-) tp — 11 

2-ti rff 






Total E-Field E j|h t (01,^ 1) & Normalized E-Field E j|h tn J0 

E4>ht( 0 1 ,4>1 ) := E4»h( 0 1 ) AFh(01 ,41 ) (e* k h cos<6l) + I'h(01 ) e‘ J khcos<6l) ) 



E*htl z := 



E<j>h 1 z - AEh 1 z - ^e"’ x ,, +rhl e 



j k-hcos(@2^ -j k h cos to 



E<t>htn( 0 1 ,<M ) ;= JMli®Lli2[ 
ma\( E4>ht I ) 



(V/m) 



(dimensionless ) 



(dimensionless ) 



(V/m) 



(V/m) 



,41) : 

(V/m) 



(V/m) 



(V/m) 
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Radiation Intensity Uh( 9 1, <j» l): 



2 

Uh(01 ,<t>l ) .=iL [(|E0ht(ei ,<M )| ) 2 + ( | E<f>ht( 0 1 ,<H )| ) 2 ] 

2t]o 



(W/solid angle) 




(W/solid angle) 



Radiated Power Pv r ^ H : 



’2 71 



Phrad -= 



Jo 



' 71 

2 



Uh(01 .*1 )■ ■sin(01)dOl d<|>l 

'0 



(W) 



Phrad = 1.059 



(W) 



Directivity Dh 0 : 

Uhmax -= Uh(0hp,O) (W/solid angle) Dho •= (dimensionless) 

Phrad 

Dho = 5.051 (dimensionless) 



Effective Isotropic Radiated Power EIRPh: 



EIRPh := Pvrad Dvo (W) EIRPh =2.133 (W) 



Radiation Resistance Rv r : 

Rhr •= — hrad ( Q) Rhr = 135.58 (Q) 

( I In | ) 2 
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Rhin = 135.58 



(Q) 



Input Resistance Rh ~j n : 
Rhin- — — (Q) 




Voltage Rejection Coeffecient fH: 

FH := — — (dimensionless) TH =0.288 

Rhin 4 Zo 

Reflection Efficiency evt- : 

eht : = ecd [ 1 - ( |th| ) 2 ] (dimensionless ) cht = 0.917 



Ga i n Gh : 

Gh = eht Dho (dimensionless) GhdB = 10 iog(Gh) 

Gh =4.633 (dimensionless) GhdB =6.659 



Polarization Loss Factor PLFh : 



Ehmax := E<t>htn(0hp,O) (V/m) 



cha : = 



0 

0 

Ehmax 



lobesh 



(dimensionless ) 



lobesh := floor 




+ 



PLFh 



( |chw cha|) 



(dimensionless ) 



PLFh = 



(dimensionless ) 



(dimensionless ) 



(dB) 



(dB) 



( dimensionless ) 
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Maximum Effective Aperture Ah ^ m : 



Ahem :=-2LL PLFh (m 2 ] 
4-71 



Ahem = 35.536 



<m 2 ) 



Maximum Effective Height hh ^ m : 



0 Rhr Ahem 
hhem .= 2 f (m) 



hhem = 7. 15 



(m) 
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E -PLANE RADIATION PATTERN: This far-field radiation pattern in the 

(x-z) plane is for the horizontal caged dipole over earth: 



Number of Dipole Conductors: 
Radius of Antenna: 

Height of Antenna Above Earth: 



N = 8 (elements ) 

a = 0 1 25 (meters) 

h = 0 25 (meters) 
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H- PLANE RADIATION PATTERN: This far-field radiation pattern in the 

(y-z) plane is for the horizontal caged dipole over earth: 



Number 


of 


Dipole Conductors: 


00 

II 

Z 


(elements ) 


Radius 


of 


Antenna : 


a =0.125 


(meters ) 


Height 


of 


Antenna Above Earth: 


h =0.25 


(meters ) 



X2(01 ) := |E$htn(01 ,0)| cos 61 + 



Y2(01 ) := | E<|)htn( 0 1 ,0)| 




0 = 0 ° 
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THE HELICAL ANTENNA (REFLECTOR OPTION) 

MATHCAD SOFTWARE-DISH JIEL.MCD 

When built to the proper specifications, the helical antenna possesses 
many qualities which make it suitable for a wide variety of 
communications applications. If the following conditions are satisfied 
the helix will exhibit a highly directional axial main lobe, low side 
lobe level, negligible mutual interference with adjacent antennas, low 
voltage standing wave ratio (VSWR) , and resistive input impedance over a 
wide frequency band: 



.8 < Cjl < 1 . 15 

n > 3 

12 < a < 3.4 

(Note: X in a subscript indicates the dimension is in wavelengths. 

Mathcad equations can not use symbolic subscripts. Therefore, the symbol 
X will immediately follow the parameter in equations (i.e., CX) in lieu 
of subscripts.) 

The helical antenna Mathcad application will compute the following 
parameters (Items with * indicate parameters that are calculated for 
both axial and peripheral feed geometries) : 

C = Circumference of Helix 

X = Wavelength 

a = Pitch Angle 

Do = Directivity 

p = Relative Phase Velocity 

Eq ^ = Electric Field Components 

¥ = Array Factor Phase Shift 

U = Radiation Intensity 

Prad = Radiated Power 

R = Antenna Input Resistance* 

T = Voltage Reflection Coefficient* 

E r = Reflection Efficiency* 
hem = Maximum Effective Height* 

G — Gain* 

EIRP = Effective Isotropic Radiated Power* 

A em = Maximum Effective Aperture* 

AR = Axial Ratio 

PLF = Polarization Loss Factor 

BW = Bandwidth 

fhigh ~ Upper Frequency Limit 
flow ” Lower Frequency Limit 
Acceptable Conductor Diameter 

E X/ y fZ ~ Electric Field Cartesian Coordinates 
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0p,(j)p = Unit Polarization Vector Coordinate Angles 

o a = Antenna Unit Polarization Vector 
r min = Mininum Distance to the Far-Field 



The following data must be input based on known or estimated data: 

D = Diameter of Helix (Center to Center) 

S = Spacing Between Turns (Center to Center) 

L = Length Along Conductor of One Turn 
n = Number of Turns 
d = Diameter of Helical Conductor 
f = Frequency of Interest 
m = Desired Mode 

i = Number of Increments in Degrees for Far Field Radiation Pattern 
I Q = Antenna Feed Current 

Z Q = Characteristic Feed Impedance 

o w = Incoming Wave Electric Field Unit Vector 



Enter input data here: 




ow ’= j 



(dimensionless ) 




0 



D 



0.1 



(meters ) 



L := 0.097 (meters) 



71 



n := 10 



(turns ) 



f := 3 10 10 (Hz) 



(dimensionless ) 




(Q) 
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S := 0.023 
d = 0.005 

x := 1 

y =i 



(meters ) 
(meters ) 

(meters ) 

(meters ) 



i =360 (degrees) 

Io = 1 (A) 

z:=1000 (meters) 



Calculate helical geometric parameters and define constants: 



c := 2.9979- 10* 


(meters/sec ) 


% - 120 * 


(Q) 


C :=n-D 


(meters ) 


X := - 

f 


(meters/cycle ) 


C =0.1 


(meters ) 


X = 9.993- 10“ 3 


(meters/cycle) 


CX := — 
X 


(dimensionless) 


js\ 

a -= atan — 

\c} 


( radians ) 


CX = 10.007 


(dimensionless ) 


a = 0.22607 


( radians ) 


LX — 
X 


(dimensionless ) 


SX :=- 
X 


(dimensionless ) 


LX = 9.70679 


(dimensionless ) 


SX =2.30161 


(dimensionless ) 


3 

O I 

00 « 

II 


(degrees ) 


k := 2 ” 
X 


(m _1 ) 


ad = 12.95276 


(degrees ) 







170 



Calculate helical antenna parameters: 



Define angular offset 6 
from helical axis: 



f, 2 n 

• s 0, — .. 2 ji 



( radians ) 



<(> := 0 , 




( radians ) 



Minimum Distance to the Far-Field 
£mi n_L 



1.61 



(meters ) 



Tj := 5-n-S 

_ _ 2 (n S) 2 
2 ' X 



(meters ) 



(meters ) 



rmin:=max(r) (meters) 



rmin = 10.58741 (meters ) 



Relative Phase Velocity p: 



P : = 



LI 



SI + m4 




(dimensionless ) 



p =2.89616 



( dimensionless ) 



Array Factor Phase Shift ^ 



y(0) := 2 n* 



|sicos(0) 




( radians ) 
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Electric Field Components E ft , E j,: 




(V/m) 



E0(0) := E(0) 



(v/m) 



E|(0) :=J -E(0) 



(V/m) 



Radiation Intensity U(0): 
U(0) •=— ( |E(0)|) 2 



(W/solid angle) 



Radiated Power P 



radi. 



Prad := 



*2 71 
0 • 



U(0) sin(0)dGdt> 



(W) 



Prad = 1.67351- 10 -3 



<W) 



Directivity D n : 

2 

Do’=12CX nSX (dimensionless ) 



Do2 = 



4 nU(0) 
Prad 



(dimensionless ) 



Do =2.7658 10 4 (dimensionless ) 



Do2 = 19.91822 (dimensionless ) 
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Axial Ratio AR : 



AR := 



LI- 




AR = 1.17586 



Effective Isotropic Radiated Power EIRP: 



EIRP := Prad-Do (W) 



EIRP2 : = Prad Do2 



EIRP =46.286 (W) 



EIRP2 =0.03333 



Polarization Loss Factor PLF : 



atan 


i-j 


[ 2 2 
'x +- y 




\ z J 



= 1.41421 10 3 (radians) 



4>p = a tan 



<|>p =0.7854 




Ex = E0(0p) cos(0p)-cos(<t>p)- E<})( Op ) • si n( 4>p ) 



Ey = E0(0p) cos(0p)-sin(<|>p) + E4>( 0p ) -cos($p) 



Ez =E0(0p)-sin(0p) - 1 



aa 



'Ex 






1 

4(|Ex|) 2 +(|Ey|) 2 +(|Ez| ) 2 



Ey 

\Ez 



oa = 



0.5- 0 5 j 
0.5 + 0. 5j 

-9.99999- 10" 4 



(dimensionless ) 



(dimensionless ) 



(W) 



(W) 



( radians } 



( radians ) 



(V/m) 



(V/m) 



(V/m) 



(dimensionless ) 



(dimensionless ) 
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Polarization Loss Factor PLF : 



,(L»I ) 2 



PLF 



(dimensionless ) 



PLF = 1 



(dimensionless ) 



Radiation Resistance R r : 

Rr:=2 - (D) Rr = 3.34702 l(f 3 (Q) 

(|lo|) 2 



Dual Parameters 



Axial Feed 
Input Resistance R: 
Ra := 140 (Q) 

Ra= 4.42874 1 0 2 (Q) 



Peripheral Feed 



Rp := 



150 

Jci 



Rp =47.41756 



(Q) 



(Q) 



Voltage Reflection Coef fecient J : 



Ra- Z, 



Ta -= 



o 



Ra+ Z n 



(dimensionless ) 



fa =0.49399 (dimensionless ) 



r . R P- Z o 

rp 

r p +z 0 

Tp =“0.51962 



(dimensionless ) 



(dimensionless ) 



Reflection Efficiency e r : 

era : = 1 - (|ra|) 2 (dimensionless ) 



crp := 1 - ( |rp|) 2 (dimensionless ) 



era =0.75597 (dimensionless ) erp =0.72999 (dimensionless) 
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Ga in G : 



Ga := era Do 


(dimensionless ) 


Gp ■= ap Do 


(dimensionless ) 


Gadb ’= 10 log(era Do) 


(dB) 


Gpdb := 10 log(ap Do) 


(dB) 


Ga =2.09088* 10 4 


(dimensionless ) 


Gp= 2.01902* 10 4 


(dimensionless ) 


Gadb =43.20328 


(dB) 


Gpdb =43.0514 


(dB) 


Ga2 '= a*a Do2 


(dimensionless ) 


Gp2 := erp Do2 


(dimensionless ) 


Gadb2 := 101og(eraDo2) 


(dB) 


Gpdb2 : = 1 0 log( ap- Do2 ) ( dB ) 


Ga2 = 1 5.05765 


(dimensionless) 


Gp2 = 14 54017 


(dimensionless) 


Gadb2 = 11.77757 


(dB) 


Gpdb2 = 11.62569 


(dB) 


Maximum Effective 


Aperture A ptn : 






A . era X 2 Do r 

Aema •= PLF 

4 71 


(m 2 ) 


Aemp ■=— PLF 

4-71 


(m 2 ) 


Aema = 0. 1 66 1 5 


(m 2 ) 


Aemp = 0. 1 6044 


(m 2 ) 


^ 1502 PLF 

4 7T 


(m 2 ) 


2 

. _ . erp-X Do2 DI r 

Aemp2 :=— PLF 

4-71 


(m 2 ) 


Aema2 = l 19657 10" 4 


(m 2 ) 


Aemp2 = 1.15545 10 4 


(m 2 ) 
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Maximum Effective Height h ^ m : 

, ^ L Aema 

hema := 2- JRr- (m) 

J % 

hema = 2 42912* 10 3 (m) 

hema2 ^Rr ^ ema j ( m ) 

4 % 

hema2 =6.51873 10~ 5 2 (m) 



hemp := 2- 

o| % 

hemp = 2.38701* 10~ 3 



hemp2:=2 |rp ^ 12171 

J % 

hemp2 =6.40574 10~ 5 



Bandwidth : 



a. u . 115 c 

fhigh - = (Hz) 

C 

fhigh =3.44759 10 9 (Hz) 

BW := fhigh - flow (Hz) 

BW = 1.04927 10 9 (Hz) 



flow := 

C 

flow = 2.39832* 10 9 



Acceptable Conductor Diameter: 

dmin := . 005-1 (m) dmax:= 05 1 



dmin =4.9965 10 5 (m) 



dma.\= 4.9965 10 4 



(m) 



(m) 



(m) 



(m) 



(Hz) 



(Hz) 



(m) 



(m) 
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HELICAL ANTENNA FAR-FIELD RADIATION PATTERN 



For the purpose of this far-field radiation pattern, the helical antenna 
axis is equivalent to the Ex = 0 grid line. The pattern is essentially 
symetric when rotated about the antenna's axis. 



Ex(8) := E 




cos(0) 



Ey(0) := E 




■sin(0) 
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PARABOLIC REFLECTOR ANALYSIS WITH HELICAL ANTENNA FEED 



The helical antenna Mathcad application (Reflector Option) will compute 
the following additional parameters (Note: a 2 or 22 will follow most 
secondary field parameters in this application, and inorder to speed up 
computation time corresponding vectors are created with subscript A & B : 

?x, y = E-Field Integral and Summation Equation Vector Components 
E ax, ay = Aperture E-Field Components 

0',p = E- Field Integral and Summation Equation Parameters 

Eq = Secondary Field Theta E-Field Component 

E^ = Secondary Field Phi E-Field Component 

U = Secondary Field Radiation Intensity 

P ra d = Secondary Field Radiated Power 

D 0 = Secondary Field Directivity 

G = Secondary Field Gain 

e ap = Aperture Efficiency 



£ s 


Spillover Efficiency 


£ t = 


Taper Efficiency 


*p 


Phase Efficiency 


£ x = 


Polarization Efficiency 


£ ohmic 


= Ohmic Efficiency 



A e = Reflector Effective Area 

h e = Reflector Effective Height 

EIRP = Reflector Effective Radiated Power 
E x , y fZ ~ Secondary Field Cartesian Components 

a a = Reflector Unit Polarization Vector 

PLF = Reflector Polarization Loss Factor 

r min = Minimum Distance to Secondary Far-Field 

0 o = Reflector Periphery Coaltitude 

0p,<(>p = Theta and Phi for Polarization Loss 



The following additional data must be inputted based on known or 
estimated data : 

E Q = Electric Field Scale Factor 

a ref = Radius of the Mouth of the Reflector 

foe = Reflector Focal Length 

ek = Blockage Efficiency 

£ sp = s P a ^ Efficiency 

tl,t2 = Radiated Power Increments for Phi and Theta 
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il = Secondary Field Increments 
N,M = Summation Increments 

r ff2 = Secondary Field Observation Distance 
x, y, z = Polarization Loss Factor Coordinates 
a w = Incoming Wave Unit Polarization Vector 



Enter additional input data here: 

foe := 5 (meters) aref := 5 

(Note: (f) should be greater 

than minimum distance to the 
far-field for the feed antenna) 



eb = .96 



esp - .95 



(dimensionless) tl =10 



t2 = 10 



( increments ) 



il =10 



N = 10 



( increments ) 



M = 10 



rfJ2 = 10 3 



(meters ) 



x = 1000 



cw2 := 




/2 



y =1000 



(dimensionless) 



z =1000 



(meters ) 



(dimensionless ) 



( increments ) 



( increments ) 



(increments ) 



(m) 



(m) 



(m) 
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Establish Integral and Summation Increments: 



*1 :=0,— ..2n 
tl 



( radians ) 



A := 1 .. N 



<f>2 :=0,— .. 2-n 
il 



( radians ) 



B := 1 M 



36, 



62 :=- 1— - icr 6 , — 2 2L_ ia 6 ..JL_ iq - 6 



il - 18 36 



36 






rp := 0,-^T. aref (m) 

t2 



*2:> := 2-it — 

B M 



rpl '•= aref — 
A N 



(m) 



022 :=|JL l.JL 

A 118/ N 36 



022 : = - — 
0 36 



Calculate reflector geometric parameters: 



aref \ 

Go := 2 atanl 1 ( radians ) 

\2foq 



r o := 1 .6-1 



00 = 0.9273 



( radians ) 



r := 10 aref 



4foc 2 +rp 

p(rp).= L. (m) 

4 foe 



r 2‘ = - 



8 aref 2 



(— ) 

\2 foe/ 



0p(rp) := 2 atanl — — ] (radians 



rmin := max(r) 



j / foe 2 * (rplj' 
P A 4 foe 



(m) 



rmin =2.0014 10 



^2 foe I 



Gpl := 2 atan ( radians 



(increments ) 



( increments ) 



( radians ) 



( radians ) 



( radians ) 



( radians ) 



(m) 



(m) 



(m) 



(m) 



(m) 
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(V/m) 



Aperture Electric Fields E A X1 E ay j_ 

Eax(rp,<t>l ) := (E0(0p(rp))-cos(<|>l ) + E4>( Bp( rp) ) sin( <J>1 )) e* k < p(rp) ” foc ' 2) 

Eay(rp, 4>1 ) = ( E0( 0p( rp ) ) sin( 4> 1 ) 4- E<K 0p( rp) ) cos(<t>l )) c ] k <p( rp) “ foc2) (V/m) 

Eaxl A B := [E9(0pl^) -(-cos(<j.l lj) -I- sinful lj ]e k ( pl A" fo 4 (V/m) 

Eayl A,B := ( E0 ( ep1 ^ 'b) + E K ep1 /') cos^ 1 ! b )) q> k ^ P ' A fOC ^ (V/m) 



Electric Field Integral Equation Vector Components P x , y j_ 
'aref 

Ax(62,*2Jl) := EaxOpJOe' k ,ps,n(92) <cos< *' >■«« *>+«*♦! >-»n<e)> -rp ^ 



(V/m) 



Px(02 > 4.2) := 



2 71 



Ax( 02 , «t»2 ,4) 1 ) dijil 



(V/m) 



Ay ( 62 , <t>2 , 1 ) 



aref 



EayCrp^D-e 1 ^m(92).(cos(0,cos(*2) + sin(*l).s,n(^2)) (v/m) 



Py(02J2):= 



2 71 



Ay( 02 , 4>2 , <(> 1 ) d4*l 



(V/m) 



Ax 1 a,b : = 



aref 

0 

' 2 71 



j k rpl sin(022/) (cos((|)ll^ cosA22n) -*-sinf^ll^) sin(<t>22^)) ... . 

Eax 1 a B -e rpl drpl (V/m) 



Pxl 



A,B 



Axl A,B d * n 



0 



Ayl 



A,B 



aref 

0 

tin 



Eayl A,B e 



k rpl sin^022^ (cos(<{>llp^ cos ($22j^ + sin sin ($22^)) 



Pyl 



A,B 



0 



Ayl AB d*" 



(V/m) 



rpl drpl (V/m) 



(V/m) 
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Secondary Field Electric Field Components E fl, 

i k e‘ j krfn 

Ee2(e2,4>2) .zL±l (Px(02,4>2)cos(i|>2) + Py(02,(|>2) sin(4>2)) 

2 -jt rff 2 



i ke‘ j krfn 

E4>2(02 , <|>2 ) -=2 ( Px( 02 , 4.2 ) - sin(4>2) + Py(02,<|i2)-cos((|>2))-cos(02) 

2 n rfl2 



• . -j k rfT2 

E <»A.B - ; J - - (^'a.B ">i* n b) * ^'a.B S '”(* 22 b)) 

2-n- ru2 

• , .j k rfTC 

E* 3 A lB r - L -^r | P,d A.B H"(* 22 b)) + fVl AiB cos(*22j ] cos(D22j 



Secondary Field Radiation Intensity U: 



U2(02 ,4*2 ) 



:= iSL [ ( | E 02( 02 ,*2 ) | ) 2 + ( | E*2( 02 ,*2 ) | ) 2 ] 
2 % 



U3 a.b~[(I e «a.bI)MI E * 3 a,bI)1 



Secondary Field Radiated Power P rP ,H : 



Prad2 := 




Prad2 =0.19318 



Secondary Field Directivity D 0 : 

Do22 l? - max(U3) 

Prad2 

Do22 = 1.04762 10 5 



(V/m) 



(V/m) 



(V/m) 



(V/m) 



(W/solid ang) 



(W/solid ang) 



(W) 



(W) 



(dimensionless ) 



(dimensionless ) 



182 



Spillover Efficiency : 



"flo 

U(0p) sin(0p) dOp 

•0 

es' = (dimensionless) 

* 71 

U(0p) sin(0p) d0p 

• 0 



es =0.19013 (dimensionless) 



Taper Efficiency gt- : 
’ Oo 



et = 



mill i ante dOp 
Prad \ 2 



0 



•2 cot 1 1 IT 



Go 



HM±f. sin( ep)d0 P 

Prad 



(dimensionless ) 



et =0.50837 



(dimensionless ) 



Phase Efficiency e^: 

- 4x6 10 ^aref 

ep:=e * (dimensionless) 



ep =0 68574 



(dimensionless ) 



Aperture Efficiency g a p £ 

ex — 98 (dimensionless) eohmic •= 98 (dimensionless) 



cap '•= es et • ep ex* eb esp ■ eohmic 



(dimensionless ) cap =0.05805 



(dimensionless ) 



Effective Isotropic Radiated Power EIRP: 



EIRP22 := Prad2 Do22 
EIRP22 =2.02379 10 4 
Ga in G : 

G22 ; = cap Do22 (dimensionless ) 

G22 =6.08192 10 3 (dB) 

Effective Aperture : 

Ae22 := eap n aref 2 

Ae22 =4.55961 



(W) 



(W) 



G22dB := 101og(G22) (dimensionless 
G22dB =37.84041 (dB) 



(m 2 ) 



(m 2 ) 



Radiation Resistance R r : 



Rr22 



Prad2 

(|lo |) 2 



(Q) 



Rr22 =0.38636 



(Q) 



Effective Height h g : 



he22 := 2 |rt22 

J % 



he22 =0.13672 



<m) 
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Polarization Loss 


Factor PLF : 




0p 2 := a , a „|ik±i' 


1 


( radians ) 


4>p2 : = atan|-^J 


0p2 =0.95532 




( radians ) 


4>p2 =0.7854 


E04 := E02(0p2, 4>p2) 


(V/m) 


E<|>4 := E«|>2(ep2,4>p2) 



Ex2 := E04cos(0p2)cos(<|>p2) - E<H s in(4>p2) 



Ey2 '= E04 cos(0p2) sin(((>p2) - E<j>4cos(<|>p2) 



Ez2 := - EG4 sin(0p2) 



'Ex2' 

oa2 — \ ■ - — Ey2 

J(|Ex2|) 2 +(|Ey2| ) 2 -»- ( |Ez2| ) 2 W, 



oa2 = 



1 0. 1 2034 + 0.48533j 
0 12034 + 0 48533j 
r 0 60556- 0.36501.) 






PLF2 : = 




2 



PLF2 =0.25003 



( radians ) 



( radians ) 



(V/m) 



(V/m) 



(V/m) 



(V/m) 



(dimensionless ) 



(dimensionless) 



(dimensionless ) 



(dimensionless ) 
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PARABOLIC REFLECTOR WITH HELICAL FEED FAR-FIELD ELEVATION PATTERNS 



E -FIELD RADIATION PATTERN: 



ll := 20 



e2 =- 1— 1 - id 6 ,— id 6 — - id 6 



36/ il 18 36 36 



X2(62) : = 



E02I02.2 



cos| 02 -t- 2 | Y2( 02 ) := 



E02|02.- 



sin| 02 + — 




- 0 . 00 .' - 0.002 -0 002 -0 001 “0 001 0 0 001 0 001 0 002 0 002 0 002 

X2(B2» 



E FIELD 
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H-FIELD RADIATION PATTERN: 



X3(62) := |E*2(02.O)| COs 




Y3i 02 ) 



|E(ii2(02.U) 





-0 005-0 002 -0 002 “0 001 “0 001 0 0 00] 0 00] 0 002 0 002 0 005 



— E FIELD 



X5< 62 ) 
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THE SPIRAL ANTENNA (REFLECTOR OPTION 
MATHCAD SOFTWARE-DISH_SPI .MCD 

Spiral antennas are a family of two or three dimensional devices that 
possess frequency independant parameters over a wide bandwidth. Spiral 
antennas are commonly used for direction finding, satellite tracking 
and missile guidance. 

The planar spiral may be of the Archemedean, log-spiral, or equiangular 
type. All three radiate two main, circularly polarized lobes 
perpendicular to the plane of the antenna. Additional gain for planar 
spirals may be achieved by placing a metal cavity on the side of the 
antenna with the unwanted lobe. The cavity may be empty or be filled 
with electromagnetic energy absorbing material. These applications 
principly examine the equiangular planar spiral and do not account for 
cavity backed effects. 

(Note: Mathcad equations cannot use symbolic subscripts. Therefore, 

symbols like /. will immediately follow the parameter in equations. ) 

The spiral antenna Mathcad applications will compute the following 
parameters for equiangular planar spirals: 

k = Wavenumber 
/. = Wavelength 

D q = Directivity 

E<)> = Electric Field Component 

U = Radiation Intensity 

Umax = Max imum Radiation Intensity 

p rad = Radiated Power 
G = Ga l n 

EIRP = Effective Isotropic Radiated Power 
A em = Maximum Effective Aperture 
BW = Bandwidth 

r min ” Minimum Distance to Far-Field 
R r = Radiation Resistance 
h em = Maximum Effective Height 
f high = Upper Operating Frequency 
flow = Lower Operating Frequency 

r n = Any point on the n ^* 1 edge of a spiral 
e ex = Equiangular Planar Spiral Expansion Ratio 
Zi = Planar Spiral Input Impedance 
T = Voltage Reflection Coefficient 
e rv = Reflection Efficiency 
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rlr = Polarization Loss Faster 

/-high = Upper Operating Wavelength 

'•low = Lower Operating Wavelength 

A = Planar Spiral Electric Field Amplitude 



The following data must be input based on known or estimated data: 



M = Mode 

N = Number of Spiral Arms 
f = Frequency of Interest 

1 = Number of Increments for Far Field Radiation Patterns 
r ff = Distance of Far-Field Calculations 
I Q = Input Current at Antenna Terminals 
r Q = Spiral Feed Point 
a = Flare Rate 

6 n+ i = Angular Arm Width of n th Spiral Arm 
<{> r = Azimuth to Compute Expansion Ratio 
R = Overall Radius 

E 0 = Source Strength Constant for Planar Spirals 

o w = Wave Unit Polarization Vector 

Z Q = Characteristic Impedance of Feed Assembly 

o a = Equiangular Planar Spiral Unit Polarization Vector 



THE PLANAR SPIRAL ANTENNAS 

Enter input data here: 



2 : 

II 


(arms ) 


f = 3.5 io'° 


(Hz) 


M = 1 


(mode ) 


lo ■- 1 


( amps ) 


Note: M max 


is N- 1 ) 






O 

m 

11 


(increments ) 


rfl •- I I0 5 


(meters 


a '= .221 


(dimensionless ) 


O 

O 

II 

O 


(m) 


70 

11 

b 


(m) 


Eo := 10 3 


(V/m) 
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(dimensionless) 



0 



8 : = 



71 



71 



3jt 

2 



( radi ans ) 



ow = 



J5 




0 



Zo := 100 



Calculate planar spiral antenna geometric parameters 
constants : 



c:=2 9979 10 8 (meters / sec) 



120 7i 



c 

- (meters / cycle) 



/. = 8 56543" 10 3 (meters / cycle) 



$ 

0 



(dimensionless) 



1 




1 C r 



»\> 



:= 4 n 10' 7 



Calculate planar spiral antenna parameters : 

Define angular offset 6 
from y-z axis: 



. 1 - io' 6 io' 6 -- icr 6 



( radians ) 



<t> :=0,- 



(Q) 



and define 



(Q) 



(Farads / m) 



(H / m) 



2 7t ( radians ) 
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Distance to Far-Field 



irr.ir.-i 



rr 0 := I 6 a 



(m) 



it j = 1 0 R (m) 



n\ 



8 R 



(m) 



rmin := ma\( rr) (m) 



rrmn =0.] 



(m) 



Wavenumber k: 




1 



(m 1 ) 



Radial Distance to n lh Spiral 
a (or - - l) 

r( n , 4>i ) •= ro e (m) 



k = 7.33552 10 2 (m" 1 ) 



r( 1.2 7i) = 4 009 1 7* 10 ' (m) 



Expansion Ratio v : 
r(n,6r+* 2 n) 

eexCn^r) (dimensionless) 

r(Mr) 

eex( 1 ,2 ;i) = 4.00917 (dimensionless) 



Bandwidth BW : 

Equiangular Spiral 
Ihigh := 4 ro (m) 



Log-Periodic Spiral 
/.high 1 := 20 ro 



Hugh := — - — 
khigh 



(Hz) 



fhighl : = — - — 
/-high 1 



Edge r: 



(m) 



(Hz) 
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thigh = 7 49475 10 ,n 


Hz 


thigh = ~ 494*5 ll) ,n 


H z , 


How := 4 R 


(m) 






flow •= ~ 
How 


(Hz) 






flow = 7 49475 1C) 9 


(HZ) 






BW .= fhigh - flow 


(Hz) 


BW1 : fhighl - flow 


(Hz) 


BW = 6 74528 IQ 10 


(HZ) 


BW 1 = 7 49475- 10 9 


(Hz) 



Electric Field E ( 0^ ) and Electric Field Amplitude A( 0 ): 
\v=U 1 (increments) 



Eo k' cos(0) ( 1 4 j a cos(0)) 3 tan 

E(0) := ; 

sin( 0 rff 




(V/m) 



A(0) =- 



cos(0)tanl-| e 



M [ — I alan( a cos( 0) ) 



sin(0) ->/] + a' cos(0)‘ 



(V/m) 



»( 0 ) 




a‘ cos(0) 2 ) 



+ atan(a cos(0) ) 



(dimensionless ) 



E<R0) := A(0) e' J 




(V/m) 
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M 




(V/m) 



Radiation Intensity U (6) : 



U(6) 



:=— (a(0) 2 ) 
2 



(W / solid angle) 




(W / solid angle) 



Umax *= max( U 1 ) 



(W / solid angle) 



Umax =2.2285-10 



(W / solid angle) 



Radiated Power P r ^n : 



n 

2 



Prad *= 4 7i 



U(0) sin(0) d0 

0 



(W) 



Prad = 7.94528- 10" 3 



(W) 



Directivity D n : 

p. _ 4 Umax , , , 

D = (dimensionless) 

0 Prad 

D q =3.52462 (dimensionless) 
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Radiation Resistance P ,1 



2 Prad 

Rr •= <Q) 



(Mr 



Rr =0 01589 



(Q) 



Incut Imoedance 



Zi 



N 30 ; 

sm| n — | 

N; 



(LM 



Zi = 1 88496* 10" <Q) 



Voltaqe Reflection Coefficient f: 



r:= 



Zi- Zo 
Zi + Zo 



( dimensionless ) 



Reflecton Efficiency c rw : 
cn := 1 - ( |f| )" (dimensionless) 



f= 0.30675 (dimensionless 1 



m =090591 ( dimensionless ) 



Gain G: 



G := cn D (dimensionless) GdB := 10 log( G) (dB) 



G = 3 19298 (dimensionless) 



GdB =504196 



(dB) 



Effective Isotropic Radiated Power (EIRP) 

Prad G 



EIRP •=■ 



(W) 



CTY 



EIRP =0 028 



(W) 



Polarization Loss Factor (PLF): 



,(uir 



PLF 



(dimensionless ) 



PLF = 1 



(dimensionless ) 
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Maximum Effective Aoertur e 



Aem ' = 



V D 0 crv PLI 
4 n 



(m 4 - ) 



Aem = 1 86416 10 



“5 



(m 2 ) 



Maximum Effective Height (h prn ) : 



, . |Rr Aem 

hem ■= J — — 2 (m) 

J % 

hem = 5 6063* 1 0 5 ( m ) 
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THE EQUIANGULAR PLANAR SPIRAL ANTENNA FAR-FIELD ELEVATION PATTERN 



For the purpose of this far-field radiation pattern, the spiral antenna 
lies parallel to the Ey = 0 grid line and is centered at the origin. 

The magitude of the electric field pattern is rotationally symmetric 
with respect to the Ex=0 grid line. The equiangular planar spiral 
antenna possesses a mirror image radiaton pattern m the -y half plane. 

E4>x(G) := |E<KG)| ccJe-t-h E<t>v(G) :* |e*( 9)| sin G+-? 




~ E FIELD 



196 



PARABOLIC REFLECTOR ANALYSIS WITH SPIRAL ANTENNA FEED 



The spiral antenna Mathcad application (Reflector Option) will compute 
the following additional parameters (Note: a 2 or 22 will follow most 
secondary field parameters in this application, and inorder to speed up 
computation time corresponding vectors are created with subscript A & B : 

P x , y = E-Field Integral and Summation Equation Vector Components 
E a x, ay = Aperture E-Field Components 

0’,p = E-Field Integral and Summation Equation Parameters 

Eg = Secondary Field Theta E-Field Component 

E^ = Secondary Field Phi E-Field Component 

U = Secondary Field Radiation Intensity 

p rad = Secondary Field Radiated Power 

D 0 = Secondary Field Directivity 

G = Secondary Field Gain 

c a p = Aperture Efficiency 

c s = Spillover Efficiency 

z t = Taper Efficiency 

Cp = Phase Efficiency 

c x = Polarization Efficiency 

c 0 hmic = Ohmic Efficiency 

A e = Reflector Effective Area 

h e = Reflector Effective Height 

EIRP = Reflector Effective Radiated Power 
E X/ y f z - Secondary Field Cartesian Components 

o a = Reflector Unit Polarization Vector 

PLF = Reflector Polarization Loss Factor 
r mm = Minimum Distance to Secondary Far-Field 

0 O = Reflector Periphery Coaltitude 

0p,4>p = Theta and Phi for Polarization Loss 



The following additional data must be inputted based on known or 
estimated data : 

a ref = Radius of the Mouth of the Reflector 
foe = Reflector Focal Length 
= Blockage Efficiency 

c sp = s Par Efficiency 

tl,t2 = Radiated Power Increments for Phi and Theta 
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ll 

N,M 

r f f 2 
x, y, 



a ref := 



foe := 



esp := 



12 := 



N := 



rfl2 : 



ow2 *: 



= Secondary Field Increments 
= Summation Increments 
= Secondary Field Observation Distance 

2 = Polarization Loss Factor Coordinates 

= Incoming Wave Unit Polarization Vector 



Enter additional input data here: 



0.32 (meters) 



0 147 (meters) cb -= 96 

.95 (dimensionless ) tl := 10 

0 (increments) il := 10 

0 (increments) M ’= 10 

: 10 3 (meters ) x := 1000 



y := 1000 

( dimensionless ) 

z := 1000 



0 



(dimension! ess ) 



(increments) 



(increments ) 



( increments ) 



(m) 



(m) 



(m) 
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Establish Integral and Summation Increments: 



*1 >0,— .. 2-ji 
tl 



( radians ) 



A := 1 N 



(increments ) 



*2 ■= Id 3 , — -+• 1(7'' 2 7i+ Id'' (radians! 



B := 1 M 



(increments ) 



62 - Id 3 ,— - Id 3 ..-- Id 

i36j il 18 36 36 



g 

d>ll D := 2 7i — (radians) 

B M 



rp :r o,^f . aref (m) 

t2 



*22:= 2n- 
^ M 



( radians ) 



rpl := aref- — 
A N 



(m) 



022 := I— I •— — — (radians) 



18/ N 36 



022 : = -- 

0 36 



( radians ) 



Calculate reflector geometric parameters: 



0o := 2 atan/- 31 ' 6 ^ 
\ 2 foe/ 


(radians) 


cT 

n 

ON 


(m) 


0o= 1.65544 


( radians ) 


r ; := 10 aref 


(m) 


. 4 foe 2 + rp 2 
W ' P) '’ 4 foe 


(m) 


8 aref 2 


(m) 


0p(rp) :=2atan| — I 
\2 foe/ 


( radians) 


rmin := max(r) 


(m) 


pl 4 foc 2 + (rp 1 A ) 2 

A 4 foe 


(m) 


rmin = 95.64028 


(m) 


h pi A \ 

0pl . := 2 atan 

A \2 foe/ 


(radians) 
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(V/m) 



Aperture Electric Fields E a YJ _ E ay _^ 

Ea\(rp,*l) :=E«K 9 p(rp)) sin(*l) V k «f’(rp»-foc 2 ) 

Eay(rp,$l ) '= E 4 < 0 p(rp))cos( 4 >] ) e* k (p<rp) “ foc 2) (V/m) 



Eaxl A.B -=E ^’( e P 1 A ) sin (' 


\ i ^ (p^A- foc ^ 
*"b) e 


(V/m) 


Ea y 1 a, B ' = 1 a) cos ( 


\ J ^ (pi a — foc 2 ) 

♦ ll B )e’ ' > 


(V/m) 



Electric Field Integral Equation Vector Components P x ^ y j_ 
“aref 



Ax( 02 ,<t> 2 , 4 >l) .= 



Eax( rp , 4» 1 ) e* k t sin(* 2 )) ^ ^ (V/m) 



0 



Px(02,4>2) •= 



Ax( 02 , 4 > 2 , 4 >] ) d<j>l 



(V/m) 



0 



aref 



Ay( 02 ,* 2 ,*l).= 

’ 2n 



Eay(rp,*l)e i k n>‘in<e2 im *»•«»< «2> + sin( ♦! >-*m «2>>. rpdrp (V/m) 



0 



Py( 02 ,$ 2 ), 



Ay( 62 ,* 2 . 41 ) d 4 >l 



(V/m) 



0 



Axl 



Pxl 



A.B 



A.B 



aref 

0 

"2-)i 



j k rpl *'"(022^ (cosAlle) cos(«2B)-l-sm(^ll B ) sm(*22 B )) ... ,,, . , 

Eaxl AB e \ n, \ \ ty \ Rr \ / \ " rpidrpi (V/m) 



AXl A.B d * n 



(V/m) 



0 



Ayl 



A.B 



aref 

0 

f 2 -n 



Eayl A ■ e j k rpl sin ( te2 ^ ( cos (*"b) C “(* 22 E() +sin((Hi^ sin^g)) ^ ^ (V/m) 



Py) 



A.B 



A y' a.b d ^ 1 1 



(V/m) 



0 
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Secondary Field Electric Field Components E q, 

i k e‘ J krfl2 

E02(02,4>2) zJ — _ ( Px( 02 . 4*2 ) cos( 4*2 ) + Pv( 02 . *2 ) sin( *2 ) ) 

2 -ti rflF2 



, -j k rff2 

E*2(02,*2) (Px(02,*2)-sin(*2)+ Py(02.*2)-cos(*2))-cos(02) 



2 n rfiQ 



E03 



A,B 



E*3 



A,B 



-j 'k' e 



-j k rff2 



2 n rffi 



, -i k rff2 

- J k e J 
2 • 7irff2 



(Pxl AB cos(«2 B )*P ) l A BS.n(42g) 

[Pxl A B (■ sin(«2 B )) * Pyl ] cos(E2j 



Secondary Field Radiation Intensity U: 



U2(02,4>2) 



:=— [ ( |E02(02 >2 )| ) 2 + ( | E*2( 02 >2 )| ) 2 ] 
2ti o 



U3 



A.B := 7~ [(| E93 A.b|) ! * (| E * 3 A,b|) 2 ] 



Secondary Field Radiated Power P rar) : 



Prad2 := 




Prad2 = 1.84163 10 



Secondary Field Directivity D n : 



Do22:=. 4nmax ( U3 > 

Prad2 



Do22 = 4.593 18 10 s 



(V/m) 



(V/m) 



(V/m) 



(V/m) 



(W/solid ang) 



(W/solid ang) 



(W) 



(W) 



(dimensionless ) 



(dimensionless ) 
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(dimensionless) 



Spillover Efficiency c g : 



es :=- 



Go 



icr 



U(Gp) sin(Gp) dGp 



-- icr 6 



U(Gp) sin(Gp) dGp 



1CT 



cs = 1.00042 



Taper Efficiency c t : 




et = 0 45135 



Phase Efficiency e^: 



-4 h-6 10 



aref 



ep :* e 



X 



ep =0.97222 



Aperture Efficiency c a p ^ 

ex •— .98 (dimensionless ) eohmic •= 98 



cap es et Ep ex eb esp eohmic (dimensionless ) cap =0.38452 



(dimensionless ) 



(dimensionless) 



(dimensionless ) 



(dimensionless ) 



(dimensionless ) 



(dimensionless ) 



(dimensionless ) 
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Effective Isotropic Radiated Power EIRP: 



EIRP22 := Prad2 Do22 
EIRP22 =8.45893- 1 0" 6 
Ga in G i 

G22'=eapDo22 (dimensionless) 
G22 = I 766)5- 10 5 { dB ) 

Effective Aperture : 

Ae22 *= cap 7i aref^ 

Ae22 =0.1237 

Radiation Resistance R r : 

Rr22 ■'= 2- Prad2 

(|lo |) 2 

Rr22 =3.68326 10 _n 

Effective Height h ^: 
he22 :=2 Rr22 

•J % 

he22 =2.19868 10 -7 



(W) 



(W) 



G22dB '= 10 log(G22) (dimensionless) 



G22dB = 52 47028 (dB) 



(m 2 ) 



(m 2 ) 



(Q) 



(Q) 



(m) 



(m) 
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Polarization Loss Factor PLF: 



0p2 := 
0p2 = 

E04 := 
Ex2 :: 
Ey2 :: 
Ez2 :: 

o a2 = 

oa2 = 

PLF2 

PLF2 



atari 




( radians ) 




( radians ) 



0.95532 (radians) 



$p2 =0 7854 (radians) 



E02(6p2,$p2) (V/m) EfJ := E4>2( 0p2, 4>p2 ) (V/m) 



E04 cos(0p2) cos($p2) - Efl-sin(^p2) (V/m) 



E94 cos( 0p2 ) • sin( $p2 ) - Eij>4 cos( $p2 ) 



(V/m) 



- E04 sin(0p2) 



(V/m) 



I 



Jl 



|Ex2|) 2 4(|Ey2| ) 2 4(|Ez2|) 2 



tx2' 

Ey2 

,Ez2 ( 



(dimensionless ) 



L 0.40672 4 0.0353 lj 
-0.40672 + 0.0353 lj 
,0.81344 - 0 07063j , 



(dimensionless ) 



• - (|crw2 oa2|) 



(dimensionless) 



= 0.16667 



(dimensionless ) 
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PARABOLIC REFLECTOR WITH SPIRAL FEED ANTENNA 



E-PANE RADIATION PATTERN: 



il :=20 




10 ' 



71 

’ll 18 



- JL- io' ; 10' 



36 



36 



Xl(02) '= 




Y 1 ( 02 ) •= 



E02 02 



,2.i) 



sin|02-«- - 




E FIELD 



Xl(02 ) 



H- PLANE RADIATION PATTERN: 



X2(02) := |E4>2(e2,0)| cos 02 + - 



Y2(02) := |E*2(02,O)| 




1.3 10 



1.18 10 



1.06*10 



94*10 



8.2 10 



Y 2(02) 7*10 



5.8*10 



4.6 10 



3.4 10 



2 . 2*10 7 



1 10 
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